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Abstract. We consider time evolution in models close to integrable points with 
hidden symmetries that generate infinitely many local conservation laws that do 
not commute with one another. The system is expected to (locally) relax to 
a thermal ensemble if integrability is broken, or to a so-called generalised Gibbs 
ensemble if unbroken. In some circumstances expectation values exhibit quasi¬ 
stationary behaviour long before their typical relaxation time. For integrability- 
breaking perturbations, these are also called pre-thermalisation plateaux, and emerge 
e.g. in the strong coupling limit of the Bose-Hubbard model. As a result of the hidden 
symmetries, quasi-stationarity appears also in integrable models, for example in the 
Ising limit of the XXZ model. We investigate a weak coupling limit, identify a time 
window in which the effects of the perturbations become signihcant and solve the time 
evolution through a mean-field mapping. As an explicit example we study the XYZ 
spin-i chain with additional perturbations that break integrability. One of the most 
intriguing results of the analysis is the appearance of persistent oscillatory behaviour. 
To unravel its origin, we study in detail a toy model: the transverse-field Ising chain 
with an additional nonlocal interaction proportional to the square of the transverse 
spin per unit length [Phys. Rev. Lett. Ill, 197203 (2013)]. Despite being nonlocal, 
this belongs to a class of models that emerge as intermediate steps of the mean-field 
mapping and shares many dynamical properties with the weakly interacting models 
under consideration. 
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The emergence of stationary behaviour in closed quantum many-body systems is perhaps 
one of the most striking features of non-equilibrium dynamics. Rather generally 
subsystems behave as they were in a “bath” and correlation functions relax to stationary 
values that can be described in a statistical fashion. 

The maturation and refinement of experimental techniques have led to the design 
of experiments ever more effective in the extraction of information on the long time 
dynamics [1-10]. A sensitive communication between theory and experiment has then 
made it possible to identify the most interesting aspects of the non-equilibrium time 
evolution in quantum many-body systems [11]. 

In particular, it was realised that integrable models behave differently from generic 
ones. The stationary properties are indeed affected by the local conservation laws, 
which in integrable models are inhnite in number. This led to the concept of generalised 
Gibbs ensemble (GGE) [12-49], which is often dehned as the mixed state with maximal 
entropy under the constraints of the local conservation laws. Non-integrable models 
with no other local conservation laws except for the Hamiltonian itself are supposed to 
“thermalise” at some effective temperature [50-56], whereas integrable models retain 
inhnite information about the initial state. 

In the same way as relaxation and thermalisation were associated with the late time 
behaviour in integrable and non-integrable models, pre-thermalisation [3,57-61] has been 
recognised as a typical feature of generic models close to integrable points. Essentially, 
at intermediate times the non-purely-elastic processes typical of non-integrable models 
are almost absent and the system behaves as if it were integrable. Despite the strenuous 
efforts to understand the process of thermalisation in the presence of pre-thermalisation 
plateaux, the picture is still far from being clear and, so far, only the earliest plateau 
has found satisfactory descriptions [59,61]. 

This state of affairs boosted the research into conserved and quasi-conserved 
operators in non-integrable models [62, 63], on the one hand, and put physicists’ 
ingenuity to the test to propose sufficiently simple models to study pre-thermalisation 
[60,61], on the other. 

One interesting proposal [60] was to break the integrability of the transverse held 
Ising chain (TFIG) by adding a highly nonlocal interaction proportional to the global 
magnetisation squared per unit length. Even though the model possesses inhnite local 
conservation laws (odd under rehection symmetry) [19,64], it was argued to behave like 
a non-integrable model in the sector of rehection symmetric states. In particular, [60] 
developed a perturbation theory that allows one to follow the time evolution of the 
ground state of a TFIG for sufficiently long times to see a pre-thermalisation plateau. 

In fact, the situation seems to be more complicated. Some techniques that will 
be developed in this paper allow us to analytically study the dynamics of a class of 
nonlocal Hamiltonians that includes the model introduced in [60]. The method we use 
is exact in the thermodynamic limit (in which the large-system limit is taken hrst): we 
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prove some conjectures of [60], but we find that the time evolution does not result in 
thermalisation. Therefore, in the thermodynamic limit, models like the one of [60] can 
not be naively related to the physics underlying pre-thermalisation. Nevertheless, we 
show that similar types of nonlocal Hamiltonians can emerge at intermediate times as 
effective descriptions of the dynamics generated by local Hamiltonians. Thus, instead of 
spoiling the interest in such models, our hndings give new motivations for their study. 

Our main goal is to investigate the time evolution of local observables under 
Hamiltonians with local interactions in the particular time windows where such effective 
descriptions can be used. There are indeed interesting cases where the expectation values 
start moving signihcantly from a plateau that could have been approximately described 
by the stationary state of the unperturbed model. We show that the crossover is driven 
by the presence of inhnitely many local conservation laws that do not commute with one 
another, which will be referred to as non-abelian integrability. We therefore identify two 
necessary requirements for a nontrivial time evolution. First, the unperturbed model 
must have a non-abelian set of local conservation laws. Second, the perturbation must 
break non-abelian integrability. 

Since the crossover appears also in the presence of perturbations preserving 
integrability, we call it “pre-relaxation” and the limit “pre-relaxation limit”. 

The pre-relaxation limit has been already considered in [65], where a typical 
crossover behaviour between two plateaux has been identihed in noninteracting models 
like the XY quantum spin chain. Ref. [65] also obtained similar results for a particular 
quench in an interacting model (XXZ spin-^ chain). However, despite an effective 
description was proposed that is supposed to capture the relaxation process for quite 
general interactions, the idea was tested only on simple cases in which the dynamics is 
essentially noninteracting. 

In this paper we start hlling this gap by investigating the pre-relaxation behaviour 
triggered off by more general (interacting) perturbations. This is a highly non trivial 
generalisation. The hrst universal picture of the time evolution of correlation functions 
after a quantum quench was delineated in conformal held theories [45,66], but most 
of the analytic results have been in fact obtained in models that can be mapped to 
free fermions or bosons [14,21,27,32-35,47,65,67-69]. For the serious complications 
introduced by the interactions, there are far less examples [70-72] in which the time 
evolution of some nontrivial observable has been worked out in interacting models. In 
the pre-relaxation limit some obstacles can be overcome. In particular we show that, 
at the leading order in the perturbation strength, the dynamics generated by (local) 
weakly interacting Hamiltonians are eguivalent to those generated by time-dependent 
(quasi-)local (mean-held) Hamiltonians, which can be solved in a self-consistent way. 
Diherently from the common situations, the mean-held mapping presented here is not 
an uncontrolled approximation, but arises naturally in the timescale investigated under 
few reasonable assumptions. The possibility to write a compact system of nonlinear 
diherential equations for the time evolution of local observables can therefore be used 
to investigate the essential aspects of the pre-relaxation limit even analytically. 
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- Section 2 is a snmmary of onr main resnlts. 

- In Section 3 we propose an effective description of the dynamics within a time 
window in which pertnrbations to a non-abelian integrable model become relevant. 

- Section 4 is devoted to identify the class of effective Hamiltonians that emerge in 
the pre-relaxation limit. We show that they can be written as polynomials of the 
local conservation laws of the nnpertnrbed model (with the correct scaling factors). 

- In Section 5 we introdnce mean-held Hamiltonians which, in the thermodynamic 
limit, generate exactly the same dynamics of the effective Hamiltonians. 

- The formalism is explicitly applied to the XYZ spin-| chain in Section 6, where pre¬ 
relaxation is investigated also in the presence of interactions that break integrability. 

- Section 7 provides a detailed analysis of the model considered in [60]. We examine 
its relaxation properties and rule out thermalisation (in the thermodynamic limit). 
Besides its intrinsic interest, the model will be useful to understand the emergence 
of oscillatory behaviour observed in the pre-relaxation limit of the XYZ model. 

- Section 8 contains our conclusions. 

- Several appendices complement the main text with the proofs of the theorems and 
additional details. 

2. Summary of the results 

We consider the time evolution of some initial state |To) with cluster decomposition 
properties]: under translation invariant Hamiltonians of the form 

H = Ho + gV, (2.1) 

where Y is a global perturbation and g is & small coupling constant. 

We focus on perturbations V that break some symmetries of Hq in such a way 
that the limit of infinite time of the expectation value of a local observable O does not 
commute with the limit of inhnitesimal g 

lim lim ^ lim . (2.2) 

g^O t^oo t^oo 

This is the typical situation in which local degrees of freedom experience a pre- 
thermalisation/pre-relaxation behaviour. Indeed, at not too large times the effect 

f We say that the state |^'o) has cluster decomposition properties if 

lim ({0i{xi)02ix2) ■ ■ ■ On{Xn)) - (Ol(xi)) {02{,X2)) ' • • (On(x„))) = 0 

min —aij |—>-oo \ / 

where the operators Oi{xi) are local (act trivially far away from the site Xi) and the expectation values 
are taken with respect to I'I'o). 
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Figure 1. Scheme of the formalism employed to investigate a pre-relaxation limit in 
models described by Hamiltonian Hq perturbed by some interacting term gV. In dark 
yellow, the section where the particular step is developed. 


of the perturbation is negligible and the expectation value has time to settle at the 
stationary value of the unperturbed Hamiltonian. On the other hand, at later times the 
perturbation can not be ignored anymore and the expectation value varies with a typical 
timescale that depends on the perturbation strength. Importantly, the amplitude of the 
variation is 0{g^), therefore the pre-relaxation behaviour can be understood in the limit 
of inhnitesimal g. 

In principle, there could be many pre-relaxation plateaux, depending on how the 
time t scales with the small parameter g. Here we focus on the limit S' <C 1 and large 
time in such a way that T = gt 0{g^). Figure 1 summarises the various steps of 
the formalism that will be developed in the next three sections to investigate the pre¬ 
relaxation limit. The analysis of explicit examples will be instead carried out in the last 
two sections. 

We now present a comprehensive summary of the main results. 

• [Section 3 and Section 4] Under some assumptions, the pre-relaxation limit of weakly 
interacting local Hamiltonians can be described by effective nonlocal Hamiltonians 
of the form 

+ + (2.3) 

where are (quasi-)local conservation laws§ of the (integrable) unperturbed 
model. In order to be a nontrivial limit, Qj ^ should not commute with one another. 

• [Section 5 and Appendix C] In the thermodynamic limit, the time evolution of a 
state with cluster decomposition properties under Hamiltonians of the form (2.3) is 
completely equivalent to the time evolution under the time-dependent mean-held 

§ It is customary to call ‘local’ a (translation invariant) conservation law with local density. 
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Hamiltonian 


liZUat) = , 

i=l j=l 

where the coefficients c^j\gt-, \&o) are obtained self-consistently as follows: 

Q(i) _ 

c^figt-,'^^) = llTT[pGGEUlA9t;'^o)^UMF{9t;^o) 


(2.4) 


(2.5) 


Here Pgge is the generalised Gibbs ensemble that emerges in the time evolution 
under the unperturbed Hamiltonian and lfMF{9t', d'o) is the time evolution operator 
of Hy[Fi9t) 

UuF{9t] ^o) = T exp drif®^(r)) ; (2.6) 

T is the time-ordering operator, which formally orders operators depending on r in 
such a way that those on the left are associated with larger (or equivalent) values 
of T. 

Despite mean-held approximations being extremely common also in the held of non¬ 
equilibrium physics [73,74], in the thermodynamic limit the mean-held mapping 
presented here is actually exact under a mild assumption. 


[Section 6] We consider Slater determinant initial states evolving under 


H = J 




1-^7 




+ 


1-7 


y^y 


-CFW. 


e'-'t+i 


+ tC/’CT, 




H— 


e ^£+2 




I e 

in the the pre-relaxation limit gt ~ 0{g^) with g 1. We identify three diherent 
behaviours: 


- Emergence of a second plateau that can be described by a GGE constructed 
with the local conservation laws of the unperturbed model. The information 
about the initial state is not encoded in a hnite number of parameters because 
the ehective Hamiltonian commutes with the unperturbed Hamiltonian and 
with inhnitely many other conservation laws in involution. If broken, one-site 
shift invariance is generally not restored. 

- Oscillatory behaviour: the expectation values of local observables keep 
oscillating with frequency proportional to g. Nevertheless, they can be 
described by a time-dependent GGE (that is not one-site shift invariant). 

- In the pre-relaxation limit the expectation values are independent of time. This 
happens whenever the initial state is one-site shift invariant, or, more generally, 
when it is an excited state of 77^p(0) (2.4) (notice that the very dehnition of 
-^MF depends on the state, so there are implicit self-consistent conditions to 
be satished). In the latter case, one-site shift invariance is generally broken. 

We stress that, rather unexpectedly, a genuine one-site shift invariant interaction 
is not always sufficient to induce the restoration of one-site shift invariance. This is 
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an indication that non-abelian integrability can snrvive interacting pertnrbations 
(at least at the lowest orders of pertnrbation theory). 

We will report explicit examples of the aforementioned behavionrs. 

• [Section 7] The eqnivalence with the mean-held description does not rely on the 
fact that the operators of (2.3) commnte with the nnpertnrbed Hamiltonian; 
relaxing snch hypothesis allows ns to constrnct more general models. Ref. [60] has 
recently proposed the model with Hamiltonian 

L \ [ ^ 

A) = - + gr,) + -[J2r,-J2 <’1 

I \ e e 

as a convenient framework for stndying pre-thermalisation/thermalisation issnes. 
Here ■ ■ ■ denotes the time average with respect to H{g, 0). In the thermodynamic 
limit L —)■ oo, we show that the time evolntion does not resnlt in thermalisation. 
The apparent conhict between onr resnlts and [60] can be traced back to the different 
order of limits. We indeed investigate the expectation valnes of local observables 
in the limit || 

(lim) lim , (2.9) 

t—Aoo L^oo 

while, for what concerns the stationary properties, [60] considered the time average 
in hnite systems 

lim-/ df . (2.10) 

T^oo T Jo 

We hnd that the stationary behavionr of local observables in the limit (2.9) is 
not characterised by a hnite nnmber of parameters. More generally, we argne 
that thermal-like behavionr can be exclnded when every linear combination of the 
operators is the Hamiltonian of an integrable model. 

• [Section 7] For Hamiltonians like (2.12) one can generally identify ‘critical regions’ 
for which the late time dynamics is nnstable nnder a small change of the parameters: 
an inhnitesimal variation can lead both to relaxation and to persistent oscillatory 
behavionr (see e.g. [73] for similar discnssions in qnantnm held theories). In 
addition, the variance of the expectation valne of some local operator in an 
arbitrarily large time window 

/1 9 /I \ 2\ 1/2 

= F” (t X * - (r X 

behaves like an ‘order parameter’ for the transition, indeed it is not analytic at the 
bonndaries of the relaxation region and vanishes inside. 

The crossover between relaxation and persistent oscillations can be illnstrated 
by means of diagrams that depict the relaxation properties as a fnnction of the 

II As a matter of fact, using the perturbative results of [60], for which the lifetime of the Ising 
quasiparticles that diagonalise Hf (7.44) scales as L~^, it is reasonable to expect that our results 
hold true even in the limit Jt ^ L“ with a < 1 and L —>■ oo. 
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Figure 2. Quench dephasing diagram of the model (2.12) in the limit of small quench 
with energy close to the ground state one. In the dark region there is persistent 
oscillatory behaviour at any time after the quench. 


Hamiltonian parameters; in particular here we focus on the limit of small quench 
(see also [21] for some clarihcations about the meaning of ‘small quench’). To 
highlight the close connection with the dephasing mechanisms [18] that allow local 
relaxation we call the diagrams ‘quench dephasing diagrams’^. We point out that 
restricting ourselves to small quenches makes it easier to interpret the results in 
terms of low lying excitations. 

We investigated the quench dephasing diagram of a simplihed version of (2.8), i.e. 

Hip, A) = - + ga’O + 1 , (2,12) 

i £ 

in the limit in which the initial state is almost the ground state of the Hamiltonian. 
In Figure 2 we can identify a ‘critical’ piecewise smooth curve at the boundaries of 
the relaxation region (the bright area). Generally, as we move towards the critical 
lines from the inside of the regions with persistent oscillations, the variance (2.11) 
of the transverse magnetisation approaches zero linearly with the distance in the 
{g,X) parameter space (c/. hgure 10). Close to the Ising critical point (^ = 1 and 
A = 0) the region of persistent oscillatory behaviour degenerates into a line ending 
at the critical point. 

Finally we show that persistent oscillatory behaviour can be related to the 

^ The terminology has not been yet standardised and in the scientific literature similar diagrams 
were sometimes called ‘dynamical phase diagrams’ (see e.g. [75]) but also ‘quench phase diagrams’ (see 
e.g. [76]) 
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3. Pre-relaxation limit 


The time evolution operator for the Hamiltonian (2.1) can be formally written as follows 
= Ui{T)e-^^°\ (3.1) 

where T = gt, 

f//(T) = exp^—i j , (3.2) 

and T1 is the anti-time-ordering operator. Uj{T) can be interpreted as the evolution 
backwards in time under the time-dependent effective Hamiltonian 

H(r) = . (3.3) 


In the pre-relaxation limit, r is hnite while g is inhnitesimal; it is therefore convenient 
to isolate the stationary (he. diagonal) contributions from V{t) 

V(r) = V + 6 V{t), (3.4) 

where V can be formally written as follows 


V = lim - 

t^OO t 


drV (r). 


(3.6) 


Le us rewrite 6V (r) in terms of its Fourier transform 6V (e) 


5V{t) = de Pf6V{e) 5V\e) = 5V{-e ). 


(3,6) 


We notice that for noninteracting Hq (which is the case we are going to consider), the 
locality of V (in the fermionic picture) implies that SV{e) is zero for |e:| > ^max 0(1). 
After some formal manipulations we hnd 

Uj{T){l - igA{T)) = I - igA{d) - i f dTUi{T)[V - ig8V{T)A{T)] (3.7) 

Jo 

with 


A{t) = i 



de eJ 


iV(e) 


(3.8) 


Eq. (3.7) makes sense as long as A{t) does. In particular, if A{t) per unit length is a 
bounded operator, in the limit —)■ 0 all the terms of (3.7) that are multiplied by g can 
be neglected, he. SV{t) is negligible. More generally, it is sufficient that the matrix 
elements of A{t) that give a relevant contribution in the expectation values of local 
observables are bounded. In [65] the operator SV{t) was worked out for a particular 
noninteracting perturbation and the previous assumption turned out to be satished. 
Besides the noninteracting case, the irrelevance of 6V{t) was also implicitly assumed 
in [61]. There, using the continuous unitary transformation (CUT) formalism [57,77], 
the authors worked out the dynamics after a quantum quench in a weakly interacting 
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model and checked the results against numerical data. For the readers familiar with 
CUT, we indeed notice that the formal simplihcation of the terms proportional to g in 
(3.7) is equivalent to CUT at 0{g^): At the lowest order of perturbation theory the 
CUT unitary transformation can be replaced by the identity; a residual dependence 
on g remains in the CUT Hamiltonian TTcut; because we are considering the limit 
of inhnite time with T = gt hnite, so Hqvt must be computed at 0{g). The latter 
is simply the time average of the Hamiltonian over its unperturbed part, namely 
Hcvt = Ho-\-gV -\-o{g). The excellent agreement of [61] with the numerical simulations 
suggests that the interaction can be replaced by its time average for rather general 
perturbations. Therefore, from now on we shall assume that the effect of 6V{t) is 
negligible in the pre-relaxation limit. 

Going back to (3.1), if 6V{t) is negligible we hnd 

e-iHt ^-iTV^-iHoT/g ^ 

where [V, Ho] = 0. We now consider the time evolution of the expectation value of some 
local operator O 

(To|e*-^*C>e-*^*|To) ^ . (3.10) 

It is well established that the stationary properties of (quasi-)local observables after 
quenches in translation invariant noninteracting models from states with cluster 
decomposition properties can be described by means of a generalised Gibbs ensemble 
(GGE) of the form 

Pgge =--, (3-11) 

where Qj are local conservation laws and Xj are real parameters determined by the initial 
state [19,34]. We also remind the reader that, in order to avoid an explicit dependence 
of the charges on the initial state, in some special cases the set of charges {Qj} could 
be non-abelian [65]. 

Let us now assume that the perturbation is sufficiently “nice” that the late time 
dynamics of under Hq can be obtained by replacing the state with the 

corresponding GGE 

{^o\e^^°T/9 ^ lim lim tr 5 [e-*^°* |To) (To| e*^“*]. (3.12) 

iS'l^oo t^oo 

This step can be easily justified for T = f if U is (quasi-)local, however in the next 
sections we’ll show that, in the presence of interactions, V belongs to a larger class of 
operators, so it is convenient to postpone the explanation of (3.12) after having clarified 
the properties of V. 

From (3.12) it follows 

lim = tr[pGGEe*^^Oe-*^^], (3.13) 

g^O 

which suggests that the pre-relaxation limit can be described by the time-dependent 
ensemble 

PtGGE(t) = e-'%GGEe‘^^S (3.14) 
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where we re-expressed the rescaled time T = gt in terms of the time. 

It is important to note that both Pgge and V commnte with the Hamiltonian. 
Conseqnently, if the two operators can be written in terms of the same set of local 
conservation laws in involntion, the time dependence disappears ptcGEit) = Pgge- In 
the next section we will show that in many cases of interest V can be approximated by 
a polynomial of the local conservation laws. Therefore, in order to see some nontrivial 
pre-relaxation behavionr, the nnpertnrbed Hamiltonian Hq mnst have a non-abelian set 
of local charges. We refer the reader to [65] for an extensive discnssion of noninteracting 
models with that property. 


4. Effective Hamiltonians 


The simplest noninteracting model that possesses local conservation laws that are not 
mntnally commnting is the XY model, whose Hamiltonian is given by 

. ( 4 - 1 ) 

i 

where cr" act like Pauli matrices on the site I and like the identity elsewhere. If the 
initial state |To) breaks one-site shift invariance, the latter symmetry is generally not 
restored in the GGE that describes local observables at infinite time after the quench. 

On the other hand, an infinitesimally small one-site shift invariant perturbation 
that breaks the non-abelian integrability of (4.1) is expected to catalyse symmetry 
restoration, which may be captured by the pre-relaxation limit. Similar issues of 
symmetry restoration have been pointed out long ago, e.g. in [78]. 

A perturbation that preserves the noninteracting character of the Hamiltonian was 
already considered in [65]. Here we investigate perturbations that have a 4-fermion 
representation in terms of the noninteracting fermions that diagonalise (4.1), namely 


V 


^l+m ^e+n2 “r+ns ^e+n4 


(4.2) 


where a“ are the Majorana fermions ({a", a^} = 25ais5in) 

= (n ^i) ^ 2/} • (4.3) 

j<e 

From the qualitative argument presented in Section 3, the relevant Hamiltonian 
in the pre-relaxation limit is determined by the time average of the perturbation. The 
calculation is not difficult but rather lengthy. However, a close inspection of the various 
contributions reveals a hidden structure that helps simplifying the computation. We 
indeed find (see Appendix B) 


Property 4.1 The time average under Hxy of a one-site shift invariant four fermion 
operator can be written as follows 


„a2 as „a4 _ 

I 




+ A 


{a} 

{n} 


(4.4) 
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where F is a linear combination of factorised terms and A is an anomalous contribution 
originated by the nontrivial solutions of the energy constraint 

s{ki) + s(k2) = ^{ks) + s{ki + /c 2 + ^ 3 ) • (4-5) 


The latter exists only in the thermodynamic limit and strongly depends on the details of 
the dispersion relation, whereas F has a structure that is almost model independent: 

1 




s=0 


ai a2 a4, 


ai as Q!2^q;4 
^713 ^712 


_l_ Q!i 0:4 ^012^0:3 

' ^ni^n4 ^712^713 5 


(4.6) 


with 


1 

L 




(4.7) 


Essentially, index s appears because the XY model (with zero magnetic held) has 
local conservation laws with momentum tt, while one-site shift invariance constrains the 
total momentum to be multiple of 27r. We also notice that the factorised part of the 
time average can be easily generalised to an arbitrary number of fermions, keeping the 
same structure of the Wick decomposition. 

To the best of our knowledge, the (quasi-)local conservation laws of the XY 
model (4.1) are noninteracting (for I 7 I 7 ^ 1) and of (4.6) seems to be a nonlocal 
conservation law that can not be (not even approximately) written as a function of 
the local charges. We then expect to become important (for local observables) 
only at times proportional to the chain length, which are far beyond the pre-relaxation 
limit. This persuaded us to conjecture that the anomalous terms are not relevant to our 
problem, which is equivalent to assume 

(To|e*^*[C>, 7l]e-*^*|To) = 0 (4.8) 


for any local observable O. In Appendix D we check the self-consistency of our 
approximation, showing that it is compatible with (4.8). We leave further investigations 
to future works. 

Proposition 4.1 suggests that in many cases of interest the effective Hamiltonian 
describing the pre-relaxation limit takes the form 




pm 


-hP ... hW + 


+ -rAF"'’ ■ ■ ■ Y™’ 

Trim-l ^ 


(4.9) 


where are (quasi-)local (he. their density is (quasi-)local, see e.g. [63, 79]) 
translation invariant operators (he. n-site shift invariant for some n G N). Indeed, 
provided that the anomalous terms in (4.4) can be disregarded, similar factorisations 
appear whenever the unperturbed Hamiltonian is noninteracting {e.g., in the model 
considered in [61]). 

Hamiltonians of the form (4.9) are therefore the perfect workbench for pre¬ 
relaxation or pre-thermalisation issues. 

We notice that the non-equilibrium dynamics generated by a subclass of 
Hamiltonians of the form (4.9) have been already worked out in [75]. The authors 
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considered ‘completely connected qnantnm models’, in which the Hamiltonian is 
symmetric nnder any permutation of the sites, and exhibited a mapping onto an effective 
classical Hamiltonian dynamics. 

We also point out that the simplest models of the form (4.9) {e.g. Curie-Weiss 
quantum Heisenberg models) have often been used as toy models to investigate the 
statistical properties in the presence of long range interactions [80]. 

The rest of the paper will be focussed on the following points: 

(i) Solution of the non-equilibrium problem for Hamiltonians of the form (4.9); 

(ii) Characterisation of the pre-relaxation limit in an interacting model, also in the 
presence of perturbations that break integrability; 

(iii) Non-equilibrium time evolution under (2.12). 

For the sake of clarity, we stress again that (ii) relies on two assumptions: 

(a) In the limit g ^ 0 with gt finite, the time evolution under H = Hq + gV can be 
split in two steps: 

1. infinite time evolution under the unperturbed Hamiltonian Wq, which is supposed 
to give rise to a generalised Gibbs ensemble |\ko) (To| —>■ Pgge; 

2. time evolution with rescaled time T = gt under the effective Hamiltonian given 
by the perturbation V averaged with respect to Hq (3.5); 

e-i^Ho+9V)t I I ^i(Ho+gV)t _ g-^^^^pGGEe*^*^ . (4.10) 

(b) The “anomalous terms” that appear in the time average of V give a negligible 
contribution (c/. Appendix D, Property 4.1 and discussion below). 

On the other hand, (i) and (iii) will be treated as ab initio problems. 

5. Solution of the non-equilibrium problem 

In this section we work out Problem (i). We are going to show that, despite the nonlocal 
appearance, operators of the form (4.9) generate a dynamics which is equivalent to that 
of a (quasi-)local time-dependent mean-field Hamiltonian. 

Here we only report some results and three useful corollaries, the details of the 
derivation can be found in Appendix C. 

For the sake of simplicity we only consider cases in which have local densities, 
however, as far as we can see, all the results can be generalised to quasi-local operators 
with tails that decay exponentially with the distance. 

In the light of (4.9), we define a class of operator £ as follows: 

Definition We say that an operator acting on a spin-| chain belongs to the class £ if 
it is written as in (4.9), namely as a finite linear combination of operators of the form 

(5.1) 

where n is finite, Hj are local translation invariant operators, and L is the chain length. 
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We consider spin chains so that the local Hilbert space is finite dimensional. This turns 
out to be a fundamental assumption for most of our results. 

One of our goals is to show that the time evolution preserves cluster decomposition 
properties, which is the key element that allows us to simplify the calculation of 
expectation values. For example we have 


Lemma 5.1 Let O E £ and |T) a state with cluster decomposition properties. The 
expectation value of O/L in IT) can be reduced to the expectation values of the local 
translation invariant operators it consists of: 


Mm 

Li —^CXD ±J J_j 



L 


( 5 . 2 ) 


3 

Using this lemma it is rather natural to relate the dynamics under (4.9) to that under 
the mean-field Hamiltonian defined as follows: 


Definition Mean-field effective Hamiltonian. Let H E £. We define the time- 
dependent mean-field Hamiltonian H^p{t) as the operator resulting from mapping any 
generic term (5.1) of the Hamiltonian (4.9) to an operator with local density, as follows 


1 




EH 

j=i 


(To|m(f)W(f)|To) 


H,, 


where U{t) is the time evolution under iL^p(f) 

U{t) = Texp(^-i driL^M^)) • 

Thus, generally must be computed in a self-consistent way. 


( 5 . 3 ) 


( 5 . 4 ) 


For example, the Hamiltonian 




X 


(5.5) 


belongs to In this trivial case Yli c’’! commutes with H, so the mean-field Hamiltonian 
is independent of time and it is given by 


Huiw = -7 E '7f is-o) E ■ 


( 6 . 6 ) 


We point out that the expectation value (per unit length) of H E £ in the state 
U{t) I To) is generally different from that of -H^p: 


1 TT 

= n {'toy'(«)^e(()i>ho) 

i=i 

n = nH (To|m(f)^t7(f)|To) (5.7) 

i=i i¥=3 1=1 


The main property that is proved in Appendix C is the exactness of the mean-field 
description in the thermodynamic limit: 






Pre-relaxation in weakly interacting models 


16 


Lemma 5.2 Let |^o) a translation invariant state with cluster decomposition 
properties and H,0 E S. Let the expectation value of O in the state that time evolves 
with H^p{t) be an analytic function oft in the strip |Im[t]| < r, with r a nonzero 

constant. In the thermodynamic limit, the time evolution with H can be replaced by the 

time evolution with the mean-field Hamiltonian: 

lim = lim o\U\t)^U {t)\^ o) ■ (5.8) 

Remark The validity of the hypothesis of analyticity on a strip can be verihed a 

posteriori. The idea is the following. The self-consistent mean-held problem can be 
generally recast into an inhnite nonlinear system of ordinary differential equations. The 
hniteness of n in (5.1) implies that the system can be written as m = F{u,t), with F a 
polynomial. If the system was hnite, the solution would have been analytic. This is not 
always the case for an inhnite system but, in practice, the numerical solution is obtained 
by introducing a cutoh parameter N that makes the system hnite. If the mean-held 
time evolution had a point of non-analyticity, the solution of the system of equations 
should display a non-trivial dependence of the mean-held parameters on the cutoh as 
iV —)■ cx). 

Corollary 5.3 Lemma 5.2 holds true in particular for local operators. 

The local equivalence with the mean-held time evolution can also be expressed in terms 
of reduced density matrices: 

Corollary 5.4 Let |To) a translation invariant state with cluster decomposition 
properties and H E S. In the thermodynamic limit, the time evolution of the reduced 
density matrix (RDM) of some spin block S is equal to the RDM in the state that time 
evolves with the mean-field Hamiltonian: 

ps{t) = l^o) (^ol = tr 5 [f/(t) l^o) (To| U\t)]. (5.9) 

The previous lemmas and corollaries are sufficient to reduce the time evolution 
under H E £ to the time evolution under a local time-dependent Hamiltonian. There is 
however another simple corollary to Lemma 5.2 that will be useful to assess whether or 
not at large times it is possible to encode the entire information about the initial state 
in a hnite number of parameters (‘thermal-like behaviour’). 

Corollary 5.5 Let H E £ and |T) a state with cluster decomposition properties. If |\h) 
is an excited state of the corresponding mean-field Hamiltonian 

= ( 5 . 10 ) 

the expectation value of local observables in |T) is independent of time. Therefore, 
|T) behaves locally as an excited state of H. 

The reverse is also true. If an excited state of H is locally equivalent to a state with 
cluster decomposition properties, then the latter is (equivalent to) an excited state of the 
corresponding mean-field Hamiltonian. 
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We are now in a position to justify (3.12), and in turn (3.13) and (3.14). 

In the limit of small g the expectation value of a local observable O reads as (3.10) 

• (5.11) 

Since Hq is local, the state I^/q) has cluster decomposition properties beyond 

some typical distance proportional toT/g (in order to be outside of the light cone). From 
Corollary 5.3 it follows that the time evolution under V is equivalent to that under the 
corresponding mean-held operator. We indeed only need JT gL {cf. (C.8)), which 
is trivially satished in the thermodynamic limit. Thus we obtain 

{^T\U'y(T)OUpT)\<tT) . ( 5 . 12 ) 


with 

Uyit) = Texp(^-i drI/MF(^)) 

and 


( 5 . 13 ) 


St) = IWo) . 


( 5 . 14 ) 


Incidentally, we notice that the time-ordering in (5.13) can not be simplihed because 
is generally written in terms of conservation laws that are not in involution with 
one another. 

For the sake of simplicity we assume that the time-dependent coupling constants 
of Cp^p are bounded. The operator U^(T)OUy{T) is then quasi-local with a typical 
range ^ proportional to T [81]. On the other hand \'^t) is the time evolution of ITq) 
at the time (cx) ^)T/g ^ T ~ which is the limit in which it is reasonable to expect 
that the state can be replaced by the corresponding generalised Gibbs ensemble (of the 
unperturbed Hamiltonian) 


(T)Oq,(T)|4-T> ~ tr|PGGEf/( (T)Oq,(r)| . ( 5 . 15 ) 

The operator V^p is obtained self-consistently by computing the expectation values of 
(quasi-)local conservation laws, which can be obtained from (5.15). Therefore, in the 
dehnition (5.3) of the mean-held Hamiltonian we can replace ITq) by Pgge 


1 




^Tj tr[pGGEC/yr)g>t/v(r)] 

i=i ^ 


(5.16) 


which is consistent with (3.13). We denote by Hmf{T) the mean-held Hamiltonian with 
the expectation values computed in the GGE. 
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In this section we investigate the (integrable) XYZ spin-| chain in the limit of small 
anisotropy in the ^ direction and also the effect of a small perturbation that breaks 
integrability. The Hamiltonian 

t i 

has the form (2.1) with Hq = Hxy (4-1) and 

’^ = 7 ^ . ( 6 . 2 ) 

£ £ 

For a fixed g ^ 0, the model is integrable for JU = h = 0, corresponding to the spin-4 
XYZ model, and for JU = 7 = 0, corresponding to the XXZ spin-4 chain; otherwise it 
is non-integrable. 

Following Sections 3 and 4, in the pre-relaxation limit (7 -C 1 with gt rNJ 0 (^°), the 
initial state can be replaced by the corresponding GGE of the unperturbed Hamiltonian 

|4^o)-tpGGE= lim lim tr 5 [e"*'^^Yt g*^^xYi] ^ 

|S|—>-oo t—^oo 

and V by the time averaged perturbation (3.5). We notice that the free Hamiltonian 
Hxy does not play any role in the pre-relaxation limit, because it commutes with pgge- 
The mapping into a mean-field problem can be decomposed in the following steps: 

- Gompute the time averaged perturbation V ; 

- Gonstruct the mean-field Hamiltonian Hmf] 

- Solve the time evolution under Hmf for any local observable. 


Some properties of the unperturbed Hamiltonian Hxy dramatically simplify the first 
step. Hxy is mapped to noninteracting fermions by a Jordan-Wigner transformation. 
Up to irrelevant (to our purposes) boundary terms, it can be written as follows 


/ aL_i \ 


H 


XY 


E 

£,n 


^2£-l ®2£-l 


^2£ 


m 


( 2 ) 

£n 


^2n- 
V 

oL-l 


V 


'' 2 n 


''2n 


(6.4) 




with a“ the Majorana fermions (4.3); PL is the block-circulant matrix 


[U] 


( 2 ) 

£n 


—e-4"-^)^-H{2)(fc), 

27r 


( 6 . 6 ) 


where the 4-by-4 matrix Pi^‘^\k) is usually called symbol (see also Appendix A) and it 
is given by 

n^P{k) = ® ; ( 6 . 6 ) 

Ek and 9k are the dispersion relation and the Bogoliubov angle, respectively 

cos k/2 -|- *7 sin k/2 


= J\ cos^ /c/2 -|- 7^ sin"^ k/2 


= 


cos^ fc/2 -|- 7^ sin^ fc/2 


(6.7) 
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Here we have chosen the two-site shift invariant representation of the Hamiltonian {i.e. 
we gathered together the fermionic degrees of freedom of pairs of adjacent sites) in order 
to be able to treat a larger class of initial states. 

In translation invariant noninteracting models almost any calcnlation can be traced 
back to operations on the symbol associated with the operator, which is the Fonrier 
transform of a block-row of the block-circnlant matrix that appears in the fermionic 
representation of the operator as in (6.4) (see also Appendix A). More generally the 
two-site representation of the symbol is a 4-by-4 Hermitian matrix, fnnction of the 
momentum and odd under simultaneous transposition and reversion of the momentum. 
A 2n-by-2n symbol completely identifies a noninteracting operator that is translation 
invariant by k sites, with k a divisor of n, by the same kind of relations that we wrote 
for the Hamiltonian [i.e. (6.4), (6.5) and (6.6)). Thus, we will often report the symbols 
instead of writing the operators explicitly. 

Coming back to the calculation of H, we find that the three constituents of the 
interaction term in (6.2) have the following fermionic representation 

\ = i > J = 1,2 , (6.8) 

e £ 

= ■ ( 6 . 9 ) 


Therefore, on the basis of our assumptions and decomposition (4.6), in the limit (7 -C 1 
with gt ~ 0{g^), we expect the local Hamiltonian (6.1) to be dynamically equivalent to 
the following nonlocal one 

H ^ H = Hxy + gV{U, h) (6.10) 


where H^y is fhe XY Hamiltonian (4.1) and the nonlocal perturbation is given by 

V(U, + BZBZ - B^BZ) + hH‘. ( 6 . 11 ) 

s=0 j=l 

The time averaged quadratic operators appearing on the right hand side of (6.11) are 
the fundamental blocks of (4.6) and read as 


1 






X(-1) 












i i 


'Z<4+j 




la^a^ 


i ^e+j ■ 


( 6 . 12 ) 












Pre-relaxation in weakly interacting models 


20 


Since we have to compute the time average of quadratic operators evolving according 
to a noninteracting Hamiltonian (c/. (3.5) and (6.4)), we can use (A.7). This allows to 
hnd the following exact result 


0{k) = lim — / dtO{k,t) = lim — 


T-S.+00 T _ ^ 

= lo{k,0) + ^^ 


T-S.+00 T 




(T^e 


0{k,d) 




a^e 


( 6 . 13 ) 


where 0{k) is the symbol of a quadratic operator and ■ ■ ■ denotes the time average. 
The symbols of the operators (6.12) read 

H:{k) =^i5s,oQ2{k)-^5s,iQsik)) 

HZ{k) = 5,^oQA{k) + 6s,iQ,{k) 

HZ{k) = 6,flQ,{k) + 5s,iQj{k) 

^"i(fc)= -^s,oQ4(fc) + 5.,iQ6(fc) 

HZi.k)= -5s,oQ3(fc) + <5«.iQ7(fc) 

^ ((1 + y) + (1 - y) cosfc) {5sflQi{k) - 5s,iQ,m 

H:i{k) = - 4 [7 + (1 - 7)(5 + cosfc)] (5,,oQ2(fc) + 5.,iQ8(fc)) 

^k 

Hll{k) = - ((1 - y) + (1 + y) cosfc) (5^,oQi(A;) + 6s,iQ^{k)) 

Hllik) = 4 (7 + s(l - 7 ) - (1 + 7)(-1)^cosA;) (5^,oQ2(fc) + ds,iQs{k)). (6.14) 

^k 

Here we expressed the results in terms of the symbols of the local charges of Hxy [65] 
Q^{k) = Xt^^\k) = Sk 
Q^{k) = Xt^^\k) = cos{k/2)ek 1 ® 

Qsik) = X^^^\k) = sin(A;) 1 ® 1 
Q^ik) = X-^°\k) = sm{k/2) 0 1 

Q^ik) = = Eu ® 

Qe(fc) = J^^^\k) = cos(A:/2) 0 a" 

Qiik) = Ji^''\k) = sin(A;) ® 

Qsik) = Ji^°\k) = sm{k/2)ek . (6.15) 


The hrst four symbols correspond to one-site shift invariant operators (the standard 
conservation laws of the quantum XY model), while the others change sign under a 
shift by one site. 

We remind the reader that from the symbol of an operator it is possible to infer its 
locality properties [65]. In particular, a smooth symbol is associated with a quasi-local 
operator. If in addition the symbol has a hnite number of nonzero Fourier components, 
as in (6.15), the associated operator is local. Equations (6.14) imply that are 
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local while are quasi-local, thus the Hamiltonian (6.10) is a member of the 

quasi-local extension of the class £ studied in Section 5. As pointed out in Section 5, 
we expect all the theorems of Section 5, in particular Corollary 5.4, to remain valid also 
for quasi-local operators. This guarantees the time evolution generated by (6.10) to be 
locally equivalent to the one generated by the following mean-field Hamiltonian 

Hmf{T) = Hxy + 2J9^((-1)‘ + 

+JgY, U’-' 

S,j 

+AAr/Zy) (6.16) 

S,j 


where {0)rp is the expectation value of the operator O in the mean-field description [cf. 
(5.15)) 


(0>T = 


Uy(T)paGEUl(T)0 


(6.17) 


To determine the time evolution generated by Hmf(T) we need to solve the self- 
consistency conditions encoded in (6.16) and (6.17). To this end, it is again convenient 
to exploit the representation in terms of symbols. Using (6.14), the symbol PiMpik^T) 
of the time-dependent mean-field Hamiltonian can be written in terms of the symbols 
{Qa(fc), a = 1,..., 8}, as follows + 


PLupik, T) = -Q,{k) + gVMF{k,T), 


(6.18) 


h 


VMF{k,T) = -^Q2{k) + ^Ca{t,ya)Qa{k) . 


a=l 


The coefficients are given by 


ci{k;yi) = 


l-|-cos/c,_(o) _(l)^ , ^ 2 ^ 


24 


■(^r+y4)+7^ 


COsfc,_(o) _(l)^ 

■(2/1 -Vi ) 


24 


~ r,^ + ^~(0) 

C2{k-,y2) = 2 —^yl 


cos k 


2U=yi^^ + 2177 


, 1 — cos k 


iV: 


f - vP 


cz{k]y^) = -U{l + '^‘^)yf^ - 17(1 - 
C4(fc; yi) = -(1 + 7^)^r - (1 - 

<=5(fc;S.) = UA7(9f+ !)«)-Vt 


cosfc ,„( 0 ) 


24 


24 


(r - 


(6.19) 


c&{k] ye) = (1 + + (1 - 7^)4^^ 

crik] m) = U{1 + 7^)y4 + ^(1 - 7^)y?^ 

/, ~ X oU — 1 9-r T cos k _(i) COS k _('o'l -('1') \ 

Csik; ys) = 27 ^—- 27 ^+ 217- 2 - ids + 2/4) 1 


( 6 . 20 ) 


■*■ From now on we set J = 1. 
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where we defined 

= 


Vo 


dp cos{£p) 
271 el 


ya{p,T) , 


( 6 . 21 ) 


and 


y^{k,T) = -Tr [Un^,{k,T)TGGE{k)Ul,^^{k,T)QM 

rT 


UHMFi.k,T) = Texp 


-i / dsVMF{.k,s) 


( 6 . 22 ) 


By taking the first derivative of (6.22) with respect to T and nsing the (closed) 
commutator algebra of Qa{k) we get 


h 


k 


ya{k,T) = ^'^fk'"'^y-f{k,T) + c^(fc;j/^)/f"'^ 2 /^(fc,T) 

/3,7=1 


(6.23) 


The nonzero structure constants that are not connected to one another by 

symmetry are given by 


fbb2 _ _ 2 fMl _ 2 

fr = fr = 2/f ^ = -24 


/r = /r = 2/r = -2(i-cosfc) 
/r = /r = 2/r = 2(i+cosfc) 


(6.24) 

The others follow from In particular Qi{k) and Q^{k) commute with all 

the other charges, so yiik) and y^ik) are conserved and the system (6.23) is reduced to 
6 first order integro-differential equations that depend on a continuous variable k. 

The solution of (6.23) entirely determines the time evolution generated by i^MF- 
Indeed, the expectation value of any local observable in the pre-relaxation limit can be 
computed using the Wick theorem with the correlation matrix 

^ ®2n-l ^ 


y 

aL-i 


V 


^2n 


^2n 


^2£-l 


''2£-l 


/ 


‘'21 


dk 


^2£ 




-i(n-£)k ^ 8yi{k,T) 


Qi{k) 


(6.25) 


^tr(Qi(A;)2) 

This also means that the reduced density matrices of subsystems are gaussian at 
any time, so the two assumptions (a) and (b) could be also reformulated as a single 
hypothesis of RDMs being gaussian. 

Equations (6.25) and (6.23) are the main results of this section: they allow us 
to compute the expectation values of local observables in the pre-relaxation limit of a 
weakly interacting model by solving a nonlinear system of differential equations, which 
is rather easy from a numerical point of view. 
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Reflection symmetry The Hamiltonian (6.1) is reflection symmetric, that is to say it is 
invariant under the transformation 


^ (^s+L-i « e {x, y, z} , 

where s is odd for reflections about a bond and even for those about a site. 
The reflection operator acts on the Majorana fermions as follows 

a? ^ 

j 

j 

Therefore the symbol "H of a one-site shift invariant operator transforms as 


(6.26) 


(6.27) 


(6.28) 


while for two-site shift invariant operators we And 




® cr^ 'H^‘^\—k) ® s odd 

0 (jy 'H^‘^\—k) 0 cr^ s even. 


(6.29) 


The symbols (6.15) of the conservation laws have the simple transformation rules 


Ql,2{k) )■ Ql,2(fc) 

Q3,4(^) 1 —Qs^flk) 

QUk) -{-^YQUk) 

Qi,^{k) ^ {-lYQj^flk). (6.30) 

Since a shift by one site is equivalent to a reflection about a bond followed by a reflection 

about a site, we recover the transformation rules pointed out below (6.15). 

If the initial state is reflection symmetric about a bond, Qj{k) = 0 for j = 3,4, 7, 8. 

Thus, the system of equations (6.23) can be reduced to 

mik^T) = -2cflk;y5)elye{k,T) + C6(fc;y6)(l + cos k)yflk,T) 
yflk,T) = -2ce{k;yG)y2{k,T) + 2(^^ -F C 2 {k;y 2 )^elye{k,T) 

y&{k,T) = 2cflk-,yY)y2{k,T) - (4 + C 2 (^;y 2 ))(l + cos k)yflk,T). (6.31) 

We numerically identified three different behaviours: 

• Stationarity. The expectation values of the observables remain equal to the initial 
values given by the unperturbed GGE (figure 3). 

• Local relaxation-. The observables relax to a different stationary value: one-site 
shift invariance is restored in some cases (figure 4) while remaining broken in others 
(figure 5). 

• Persistent oscillations: The amplitude of the oscillations of the expectation values 
of the observables does not approach zero (figure 6). 
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gt 


Figure 3. The time evolution of (red dashed) and (cf erf ) (blue) after a quench 

from the state |MG) (6.43) and Hamiltonian H (6.1) with 7 = 2, h = 0, and U = 5. 
The correlators are stationary. We find stationary behaviour whenever the initial state 
is reflection symmetric, y 2 {k) = ye{k) = 0 (c/. ( 6 . 22 )), and h = 0. 


We point out that, even when there is relaxation (at some intermediate times with 
Jt 3> g~^), the stationary state is not thermal, being the local conservation laws of 
Hxy with symbol proportional to Qi{k) and Q^{k) (namely the charges that preserve 
non-abelian integrability) conserved in the pre-relaxation limit. 


6.1. Perturbations preserving integrability 


In this section we consider the case f/ = 7 = 0, in which H (6.1) is the Hamiltonian of 
the XXZ spin-| chain. Because of the U{1) symmetry of rotations around 2 ;, there are 
many simplihcations and the system of equations (6.31) can be rewritten as follows: 


kypT) = y'^(T)y'^(T) - vf\T)vt\T) 

\yt\T) = -(2yt\T) + y'r'-\T) + y7‘l(r))!/?(r) + (ft + 4y!,“>(T))!,l”l(r) 
\yt\T) = + !/7‘'(T))t/f'(T) - (ft + 4y!,”'(T))!,'"l(r), (6.32) 


where we dehned 


yt\T) 



dp cos (np) 
27r 1 -b cos p 


1/2 (p, T ), 



dp 


cos(np)?/5,6(p,T). 


(6.33) 


We notice that, despite the denominator, are expectation values of local operators, 
as well as y^l- Since \ commutes with the Hamiltonian, the dependence 

on h is simple and, in particular, the expectation value of the one-site shift invariant 
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gt 


Figure 4. The time evolution of (erferf) dashed) and (erferf) (blue) after a 

quench from the state |MG) (6.43) and Hamiltonian H (6.1) with 7 = 2, h = 1 
and U = —2. The correlators rapidly relax to the same stationary value, restoring 
translation invariance. We verified relaxation up to gt = 1000. 



gt 


Figure 5. The time evolution of (erf erf) (red dashed) and (erf erf) (blue) after a quench 
from the state |MG) (6.43) and Hamiltonian H (6.1) with 7 = 2, h = 2 and U = 2. 
The correlators rapidly relax to different stationary values. We verified relaxation up 
to gt = 1000 . 
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Figure 6 . The time evolution of (crfa^) (red dashed) and {a^crf) (blue) after a quench 
from the state |MG) (6.43) and Hamiltonian H (6.1) with 7 = 4, ft, = 1 and U = —2. 
The correlators exhibit persistent oscillations on the time window explored. Inset: the 
amplitude of the oscillations is still unabated at gt = 1000 . 


conservation laws is independent of the magnetic field. This means that the fnnctions 
t/ 2 ”^ are independent of h. 

It is usefnl to rewrite the system for n = 0. We hnd 
yf{T) = 0 

4°'(T) = 2{h + 2yf{T) - 2yf{T))yf{T) 

yf{T) = -2{h + 2yf{T) - 2yf{T))yf{T). (6.34) 

Inspecting the system we conclude that and ( 1 / 5 '')^ + are conserved. Moreover, 

the system (6.34) can be directly solved, it yields 

( 5 ^) ^ ^ 

2 


y?{T) = 


2L 


yt\T) = 2/?'(0) cosf / dr {2h + 4m(r))) + t/|°^(0) 


sm 


dr {2h + 4m(r)) 


yf^{T) = yf^{t]) cosi^j dr (2h + 4m(r)) j — t/|’^(0) sin^ / dr (2h + 4m(r)) j . (6.35) 


Here we defined 

{rji\ _ ^ h] 

^{T) = — - yl^\T) = s - — 


= T ) + r 




(6.36) 


If both 2/5 (T) and t/f* (T) are zero, (6.32) has the solution 
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yP{T) = 2/^(0) cos(2(/i + 2s^)T) + 2/["](0) sin(2(/i + 2s^)T) 

yf^ (T) = y'f (0) cos(2(/i + 2s^)T) - yf^ (0) sin(2(/i + 2s^)T). (6.37) 

For h 7 ^ —2s^, local observables keep oscillating in time, otherwise, on the pre-relaxation 
timescale, the expectation values of local observables do not move from the values 
reached at times 1 Jt g~^. 

More generally + ( 2 / 6 °^)^ 7^ O)’ from (6.35) and (6.36) we immediately infer 

that relaxation is possible only if 

3 1imm(r) = -^, (6.38) 

T^oo Z 


3 lim 

T^oo 



< (X) . 


(6.39) 


We see that m(T) could be interpreted as a sort of ‘induced magnetisation’ that h must 
compete with. 

The trivial dependence on h is manifest choosing the variables 


Yn = v'l'm 

4,. = !/!,”’ (T)j/f (T) + yf' (T)yf (T ), (6.40) 


which satisfy the following system of equations independent of h: 

m(r)^|-yi(T) 

W(T) = - 4 $o( 2 F;(T) + W+i(T) + y;_i(T)) -l- 8 ( 5 ^ - m(T))<F„(T) 

$„(T) = -2(s^ - m(T))W(T). (6.41) 


Here we omitted the time dependence in the conserved quantity <Fo. Since $0 7 ^ 0 by 
assumption, the original variables are obtained from the inverse transformation 


, [n]_. _ 2|/f(T)$„(T)+i/[°l(T)W(T) 
^ > 2 <Fo 


yT\T) = 


2yf\T)^^{T)-yf\T)Yr,iT) 


2<Fr 


(6.42) 


and (6.35). Performing a qualitative analysis of the system (6.41) we conclude that 
condition (6.38) and m{T) + h/2 approaching 0 faster than 1/T imply relaxation of 
local degrees of freedom. Therefore the variance (2.11) of m{T) is what in Section 2 
we called an ‘order parameter’ for the transition between relaxation and oscillatory 
behaviour. 

As we will show in Section 7, some aspects of the solutions of nonlinear systems 
like (6.41) can be worked out analytically. In the present context this would involve 
the study of quantum quenches from rather artihcial initial states. Therefore we prefer 
to leave the entire discussion to Section 7, where we will obtain a system of equations 
extremely similar to (6.41), with the advantage that the qualitative analysis can be 
carried out for more conventional initial states. 
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In order to gain some insights into the time evolution under the Hamiltonian (6.1) in 
the non-integrable case we focus on quantum quenches starting from the dimer product 
state 

iti) - w) „ „ itw - in> 


|MG)= 2^-2. 

which is the ground state of the Majumdar-Ghosh Hamiltonian 

L 


TT ^ 

^ CTe ■ CTi+X + -CTl ■ <7^+2 • 


(6.43) 


(6.44) 


i=\ 


Despite the model being interacting, (6.43) is a two-site shift invariant Slater 
determinant, whose correlation matrix has the following symbol 

V^G{k) = a^®aP (6.45) 

The initial conditions for {ya{k)} (6.22) are determined by the GGE correlation matrix 
for (7 = 0. They can be obtained by expanding rMG(^) in fhe base of the symbols (6.15) 
of the conserved charges of Hxy (the remaining space is zeroed by the time evolution, 
as 1 ^ Jt, cf. section 5.1) 




i=l 


tr[(Qi(fc))^ 


1 -|- cos k 


-Qi{k) - 7: 


1 — cos k 


We find 

’ 1-f-cos fc -I- 7^(1 — cos fc) ^ ' 1-|-cos/c-|- 7^(1 — cos/c) 

The only nonzero initial conditions are given by {cf. (6.22)) 

1 -|- cos k 


(6.46) 


Qs{k). (6.47) 


2/i(fc,0) = 

2/5(fc,0) = -7 


4 

1 — cos k 


(6.48) 


The initial state is reflection symmetric about a bond, so we can use the reduced system 
(6.31). Since y^ is the only nonzero initial condition that appears in (6.31), for 7 = 0 
(see previous section) the solution of the system of equations is independent of time, 
namely the pre-relaxation limit is trivial. 

It is easy to see that also for h = 0 system (6.31)(6.48) has a stationary solution. 
We therefore assume 7 , h 7 ^ 0. Since Cj{k;yj) are linear homogeneous functions of yj, 
the magnetic held h enters into the equations essentially as a scale factor. We rescale 
the variables as follows 

7 2yj 2cj 

Yh = - ^^ = -- 


e = 


T = 2hT = 2hgt 


(6.49) 
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drZ2{k,T) = - eel'y5{k,T)ze{k,T) + ecos'^ -'ye{k,T)z5{k,r) 

drZ5{k,T) = - e'yG{k,T)z2{k,T) + (1 + eel'y2{k,T))z6{k,T) 

COS^ — 

d-rZeik^r) = e-fG{k,T)z2{k,T) - (1 + eel'^2{k,T)) — ^^z^ik^r ), 


with the initial conditions 

Z2{k, 0 ) = z^ik, 0 ) = 0 z^ik, 0 ) = 


et 


. 9 k 
■ sin^ -. 
2 


(6.50) 


( 6 . 51 ) 


For generic e the system of equations is not exactly solvable, but the limit of small e 
allows a linear approximation. For not too large rescaled times (we’ll come back to 
this point later) the terms that are multiplied by e in the last two equations can be 
neglected, while the functions that appear on the right hand side of the first equation 
can be computed at 0 (e°). 

For zg and zq we obtain the simple solution 


z^ik^r) 


k 


sm - cos ( 


/cos 


-r 


V Sk 


Ze{k,T] 


o k cos ^ 

^ sm —-^ sm 


/cos 


2 ek 


-r 




(6.52) 


while Z2 is a slightly more complicated function that involves integrals over the 
momentum of 2:5^6) namely 

sin A; sin I r dp ^ > 

' 2^9aik,p)f{- 


Z2{k,T) ^ e 
Here we defined 


27r 


,cos I cos 2 
- + CT- hr) 


7 =± 




sin^i 


,2 k 


P P 


,2 k 


Sk 


P P 


( 6 . 53 ) 


cos"' f cos"' § — 7 " sm" ^ sm" § k p 

- - - - -^- - + a cos - cos - 

2 2J 


Sk£p 




1 — cos(a;r) 


X 


For small e and given 7, Z2{k,T) relaxes to the stationary value 


Z2{k, 00) = 


e sin^ k 
872 


1 — 37^ + (2 + 47^) cos /c + (1 — 7^) cos 2 k 


(6.54) 

(6.55) 

(6.56) 


Let us now estimate the time window in which the linear approximation is applicable. 
From (6.52) it follows that 75 and 70 decay to zero as r“ 2 . Instead, since Z2 approaches 
a nonzero stationary value, 72 is of the same order of Z2- This means that, as the 
time increases, the first term on the right hand side of the last two equations of (6.50) 
becomes more and more negligible with respect to the other term multiplied by e. By 
neglecting the former we obtain essentially the same solution (6.52) as before, with the 
replacement 


r -)■ r + e 


ds72(A;, s) = r (1 + e el'j2{k, 00)) + ... . 


(6.57) 
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<Q5> 
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gt 


Figure 7. The time evolution of =7 / where Qs is the conserved 

charge of Hxy corresponding to the symbol Q^lk) (6.15) for a time evolution starting 
from the state \MG) (6.43). The parameters of the Hamiltonian (6.1) are 7 = 0.2, 
h = 3.5 and U = —1, hence e « 0.029 {cf. (6.49)) and 7 fulfils the consistency condition 
(6.58) of the linearisation procedure. The analytical prediction of (6.52) (red dashed 
line) is in excellent agreement with the numerical data (blue line). 


Being 72 ~ 0(e), after a rescaled time t ~ the correction to becomes comparable 
with the function itself. Assuming that the relevant part of the time evolution occurs 
within this time scale, the linear approximation is justihed only if 1^21 <C 1 (and e <C 1 ). 
For 7 < 1/2 we hnd \z 2 {k,oo)\ < so we obtain the consistency condition 

^«7«2h. (6.58) 


Figures 7 and 8 report a comparison between the solution of the linearised problem and 
the full numerical solution of system (6.31) for a set of parameters fulhlling (6.58). 

From the expressions (6.52) of Z 5 {k, r), ZQ^k, r) and (6.56) of Z 2 {k, r), we can directly 
compute the time evolution of the expectation value of any local observable in the pre¬ 
relaxation limit. Indeed Corollary 5.4 allows us to apply the Wick theorem at any time 
(in the limit under examination) and the correlation matrix is given by (6.25). 

Any integral involving z^ and Zq approaches zero and Z 2 becomes independent of 
time even if not integrated. Therefore, in the limit (6.58) and for large times the 
expectation value of any local observable relaxes to a stationary value that can be 
described by the correlation matrix with symbol 


'yzi{k,oo) 'yz2{k,oo) 

hm F(fc;r) = - 2 -^’ 


r^oo(*) ' ' ' e% ' el cos'^{k /2) ^^' (6.59) 

where the inhnite time limit r —)■ 00 ^*^ must be understood within the limits of validity 
of the linear approximation. 
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Figure 8. The time evolution of {Q2)t = if where Q 2 is the conserved 

charge of Hxy corresponding to the symbol Q 2 {k) (6.15) for a time evolution starting 
from the state \MG) (6.43). The parameters of the Hamiltonian (6.1) are 7 = 0.2, 
h = 3.5 and U = —1, hence e « 0.029 (c/. (6.49)) and 7 fulfils the consistency 
condition (6.58) of the linearisation procedure. The analytical prediction of (6.52) 
(red dashed line) is in fairly good agreement with the numerical data (blue line). The 
stationary value produced by the solution of the linearised problem (black dotted line) 
is (Q 2 ) = ^\~ 28 l ^ (c/- (6.56)). 


We point out that one-site shift invariance is restored, indeed the only contributions 
to the correlation matrix at inhnite times arise from Qi{k) and Q 2 {k), which are symbols 
of one-site shift invariant operators. The manifestly one-site shift invariant expression 
of the correlation matrix in the limit (6.59) reads 

with 

2h(l -b cos 2k) -b sin^ k cos k{l — 87 ^ -b (2 -b dy^) cos 2/c + (1 — 7 ^) cos 4k) 

2h{l -b cos 2k + 7^(1 — cos 2k)) 

X {cos k — 1 sink a^). (6.61) 


6.3. Remarks on the late time dynamics 

We notice that, despite one-site shift invariance being restored in (6.59), the asymptotic 
value is not given by the average over a shift of the expectation value of the operator 
in the GGE of the unperturbed model. Indeed the one-site shift average of (6.47) is 
proportional to Qi{k) {cf. (6.30)) but the symbol of the large time correlation matrix 
(6.59) has also a term proportional to Q 2 {k). Gonsequently, the shift-averaged stationary 
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values can not be recovered from those in the limit of small perturbation 0. For 
example we have 



(6.62) 



(6.63) 


lim 

2hgt^oo^*) 


where we highlighted that there are 0{g) corrections. Besides this particular quench 
in a non-integrable model, similar issues arise also in the integrable case, where it is 
generally believed that at infinite time after the quench the expectation values can be 
computed in a GGE constructed with the (quasi-)local conservation laws of the model. 
In the rest of the section we propose a reasoning that relates this kind of deviations to 
possible atypical properties of the model. 

We notice that at times 1 f shift-symmetrised expectation values of 

quadratic operators, e.g. (6.62), can be obtained from the GGE of the unperturbed 
model constructed with only the local translation invariant conservation laws. This 
is because the limit 1 ^ f -C - with —)■ 0 can be described by the GGE of the 
unperturbed model, but the conservation laws that break translation invariance are odd 
under a shift by one site, so they can in fact be neglected (this equivalence breaks down 
for operators that consist of the product of more than two Jordan-Wigner fermions). 

For nonzero g non-abelian integrability is supposed to break down and the relevant 
charges are generally assumed to be one-site shift invariant and in involution with one 
another. We now speculate about the stationary state in the limit of small g if the 
perturbation does not break integrability. Let the (quasi-)local conservation laws of 
the interacting (integrable) model be in a smooth one-to-one correspondence with the 
local one-site shift invariant conservation laws for g = In the limit of small g, the 
stationary state should locally approach the GGE constructed with the local one-site 
shift invariant conservation laws of the unperturbed model. The expectation values of 
shift-symmetrised quadratic operators at times 1 f ^ are compatible with such 
a one-site shift invariant GGE. However, discrepancies like that between (6.62) and 

(6.63) show that at larger times there is a time window in which the expectation values 
approach a different value. Our assumption of regularity of the conservation laws as a 
function of g requires that at even larger times the expectation values should eventually 
relax to the same values they had in the earliest plateau. This is clearly possible, but 
an inhnite number of operators displaying a similar behaviour is rather surprising. This 
makes us wonder whether the hypothesis of regularity could break down, that is to say 
there are (quasi-)local conservation laws for g ^ Q that become nonlocal when g = 0 
{e.g. their typical range could be singular as —)■ 0) or, vice versa, some one-site shift 
invariant conservation laws of the unperturbed Hamiltonian do not have analogues at 
nonzero g. 

From this point of view, discrepancies like that between (6.62) and (6.63) could be 
indications that in the XYZ model there might be quasi-local conservation laws that do 
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not behave well in the limit g ^ 0. This scenario becomes even more plausible if one 
takes into account the issues in the construction of the GGE in XXZ spin-^ chains that 
were recently pointed out [42-44], 

6.4- Summary 

In this section we considered the pre-relaxation limit in an XY spin-| chain perturbed 
by interacting operators. 

As an example of pre-relaxation in integrable models we investigated the XXZ 
spin-^ chain. The model has U{1) symmetry, manifested by the conservation of the 
spin in the 2 ;-direction. Gonsequently, the external magnetic held (along z) generally 
produces oscillatory behaviour in local observables. In fact we showed that there is a 
specihc (generally nonzero) value of the magnetic held for which the time evolution in 
the pre-relaxation limit may end up in a second plateau. 

For non-integrable perturbations we exhibited examples of the typical time 
evolution of the expectation values of local observables in the pre-relaxation limit. 
We also described a linearisation scheme that allowed us to predict the time evolution 
of the dimer product state (6.43) when the Hamiltonian parameters satisfy particular 
conditions. In that limit we found local relaxation to a one-site shift invariant state. 
In order to characterise the crossover between persistent oscillatory behaviour and 
relaxation, one should go beyond that linearisation scheme. 

In the next section we will consider the model (7.1), which has several dynamical 
aspects in common with the Hamiltonian (6.10), especially in the integrable case 
7 = t/ = 0 considered in Section 6.1. Specihcally, we will present a method that 
resembles the linearisation considered in this section but that allows us to extract some 
information about the ‘quench dephasing diagram’ of the model. 

7. An exactly solvable model 

The results of Section 4 are a compelling motivation for the study of nonlocal 
Hamiltonians of the form (4.9). We now go beyond that rigid derivation: we skip 
the formal steps that relate (4.9) to a pre-relaxation limit and start off directly with a 
Hamiltonian of the form (4.9). We then query whether such models with (apparently) 
non-integrable long-range interactions could display thermal-like behaviour at late times 
after a quench. Specihcally, we consider the Hamiltonian 

h{3, A) = - + acl) + ^ (E <^l)7 (7.1) 

e e 

This has been recently proposed as a convenient model to investigate pre-thermalisation 
issues [60]. In fact, in order to recover some temporal cluster decomposition properties, 
the authors of [60] considered a slightly modihed version of (7.1), where, in the term 
proportional to A, replaced by 'Y^i~ latter being its time 

average for A = 0. The time average is a simple quadratic quasi-local operator (c/. the 
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first equation of (6.14)), so our formalism could be readily applied. However, for the 
sake of simplicity, we consider (4.9) and show later that a redehnition of g is sufficient 
to recover the results shown in [60]. 

From Corollary 5.4, in the thermodynamic limit L —)• oo the time evolution under 
(7.1) is locally equivalent to the time evolution under the time-dependent mean-held 
Hamiltonian 

+ m<rl ). (7.2) 

I 


with hit) the solution of the self-consistent equation 

hit) = ^ - 2A (4'o|T'^exp(^z^ dri7M^(r) Texp(^-i drif^U^) l^o) • (7.3) 

Here we used translation invariance to replace by the local operator (which 

removes the nasty dependence on the chain length L). For Slater determinant initial 
states the expectation value can be conveniently written in terms of the symbol F(/c) of 
the initial correlation matrix (see Appendix A). In particular we hnd 

/ 77 AU 

(7.4) 

with 


Vkit) =tr UnMAk,t)T{k)Ul^^{k,t) 




(7.5) 


f7«Mr(*:. *) = Texp 2i / dTtr''(/i(r) - e*'”’”) 
One can easily verify that is the solution of 


ylit) = 4(h(t) — cos k)(j)k{t) + 16(h(t) cos k — 

= -4h(t)|/l(t), 


with 


l)!/s(*) 


(7.6) 


(7.7) 


0fe(O 


—4tr 




(7.8) 


Equations (7.4) and (7.7) could be solved by discretising the momenta k; however 
working in real space is more transparent. We therefore introduce the (real space) 
Fourier coefficients 


Un 



dk 


cos{nk)ykit) 





cos{nk)(j)k{t ), 


(7.9) 


which are the expectation values of local operators with range n and n -|- 1 respectively. 
The system of equations (7.7) can then be recast as follows 


h = 9-Xyo 

y^ = Ah(j)Q - 401 -f 16h^i - 16^o 

= Ah^n - 2(0n+i + 0n-i) + 8h(y„+i y„_i) - 16^„ {u > 0) 

= -^hy'n (n > 0). 


(7.10) 
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The similarity with the system of equations (6.41) for the pre-relaxation limit in XXZ 
spin chains is remarkable, although the meaning of the variables is different. 

For sufficiently large n (in order to be outside of the (deformed) light cone), both 
ijn and are small (in the noncritical case they are exponentially small in n). Thus, 
the error originated from truncating the system of equations to n < N decays very fast 
to zero with N and (7.10) can be conveniently reduced to a finite system of differential 
equations. This is the regularisation that we used in our numerical investigations*. 

The system of equations has at least one integral of motion, namely the energy per 
unit length £ = {4/o|i7|\I/o) jL. This can be written as follows 



(7.11) 


In addition, (7.1) has infinite noninteracting conservation laws that are odd under 
reflection symmetry [19,64]. The Dzyaloshinskii-Moriya interaction 



(7.12) 


is one of them. For generic initial states, this is sufficient to rule out thermalisation. 
However, we embrace the point of view of [60] and wonder whether at infinite time after 
a quench from a reflection symmetric state some form of thermalisation arises. 

7.1. To relax or not to relax 

At late times the mean-field Hamiltonian (7.2) can result in two distinct behaviours: 
either the dynamics is governed by a (time independent) TFIC Hamiltonian (viz. 
31imt_>.oo h{t)), or there is no relaxation (viz. Hlimj^oo h{t)). 


In the following we will focus on quantum quenches from the ground state of the 
TFIC Hamiltonian 




(7.13) 


e 


A numerical analysis suggests that for generic values oi go, g and A, the system of 
equation (7.7) does not always describe a relaxation process. We indeed found cases in 
which h{t) oscillates with almost constant amplitude (see Figure 9). 

If there are (finite) regions of the parameter space associated with relaxation and 
regions that are instead characterised by persistent oscillations, some quantities shall 
behave non-analytically at the boundaries of the regions. For example, the relaxation 
parameter 



(7.14) 


* From a numerical point of view, this allows us to avoid integrating at each time step of the Runge- 
Kutta algorithm used for the resolution of (7.10), at the cost of complicating the initial conditions. This 
is convenient because generally the number of time steps is much larger than N. 
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Figure 9. The time-dependent magnetic field h{t) of the mean-field Hamiltonian (7.2) 
(essentially, the magnetisation along z, cf. (7.3)) after a quench from the ground state 
of the TFIC (7.13), with go = 1.5, and Hamiltonian H{g,X) (7.1), with g = 0.5 and 
A = 0.6. The system does not seem to relax, indeed also at very large times (inset) 
there are (rather regular) persistent oscillations. 


approaches zero as T —)■ cx) in the regions of (local) relaxation. On the other 
hand, if there are (sufficiently regular) persistent oscillations, AIit remains nonzero 
for arbitrarily large T. 

We numerically analysed the region of the parameter space in which there is no 
relaxation (at least apparently). This generally happens for sufficiently large A (see 
Figures 2 and 9). In the vicinity of (a numerical estimation of) Ac, AIit is nicely htted 
by a line {cf. Figure 10) 

AhT = a(A-Ac) + 0((A-Ac)^) A>Ac. (7.15) 

Since AIit is positive, linear behaviour is indicative of discontinuous derivative at A = Ac- 
Indeed, for A < Ac, Ah^ is compatible with zero (see e.g. the three solutions with A < Ac 
in Figure 10). 


7.1.1. Small guench. In the case go = h{0) the solution of (7.10) is independent of 
time (as a consequence of Corollary 5.5). This can be turned into a condition on the 
parameter g of the Hamiltonian as follows 


9 = 9{9oA) , 
where we dehned 


dk 


go — cos k 


(7.16) 


fi'(fi '05 — fi'o + 2A 


-TT 2vr ^1 + 5(2 - 2go cosk 


(7.17) 
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Figure 10. The relaxation parameter Ahx (7.14) for quenches from the ground state 
of the TFIC (7.13), with go = 1.5, and Hamiltonian H{g,X) (7.1), with g = 0.5, as a 
function of the rescaled time Ac+ 7^(A —Ac) for various values of A. The maximal time 
considered for the time average is Tm = 1920; at Tm the abscissa is exactly equal to 
A, which can therefore be identified with the abscissa of the open squares. The critical 
value Ac « 0.5257 (grey vertical line, where all curves collapse) has been estimated by 
a parabolic fit (black solid line) of Ahx (open squares) for large T as a function of A. 
The dashed line is the linear term of the fit. 


The initial state {ipgo) corresponding to this quench is then an effective eigenstate of 
H{g,X) (7.1). In addition, in Appendix E we show that one of the solutions of (7.16) 
corresponds to the state iV'go) with minimal energy among Slater determinants. As a 
matter of fact, the numerical analysis indicates that \'ipgo) is the true ground state of 
(7.1). We can therefore use I'lpgo) as a reference state to dehne the limit of small guench. 

In this limit, the transition relaxation/no-relaxation can be understood more clearly. 
Indeed, choosing the parameters such that (7.16) is approximately satished, both 
IVnit) — ^n(0)| and \(pn{t) — 0n(O)| tum out to be small at any time. We can therefore 
linearise the system of equations (7.10) isolating a time independent contribution from 
Vn snd (pn- 

h = go+ {g- g{go, A) - X6yo) 
yn yn T dyn 

0n 0n T dcpYi 1 

where variables with a bar on top are expectation values calculated in iV'go)- AVe then 
obtain 


n 1 


“ -16(1 + al)iyn + 169(](i59n+i + hn-O “ >'vjya + (g - g)v, 


(7.19) 
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where 


'^n 4(0,2 -|- 2yYi—\ 2^2j_|_]^) . 

(7.20) 

The system of equations (7.19) can be readily solved, 
ground states of TFIC Hamiltonians y'niO) = 0, we hnd 

Since for quenches from the 

Syn{t) ^ {g - g)'^kj 1 - cos{y/a~t) Wnj 

(7.21) 


j 


where aj and Wnj are the eigenvalues and the components of the (right) eigenvectors 
(at hxed j) of 

A^n 16(1 9o')^^n 165^0*^|£—n|,l A (7.22) 

and kj are given by 

k = {AW)-^^ {[k]j = k,, [W]nj=w^,,, [^j = vj). (7.23) 

For A = 0, A can be diagonalised in momentum space; the eigenvalues are given by 

ttj = 16(1 + gl- 2go cos kj ), (7.24) 

and, in the limit N ^ oo {N is our regularisation parameter), the spectrum becomes 
continuous (the density of kj is uniform in (0, tt)) and the eigenvectors unbounded. It is 
not difficult to show that the rank-1 perturbation Xv^Sno does not change the continuous 
part of the spectrum, which is still described by (7.24) and, as —>■ oo, the density of 
kj remains uniform. 

For sufficiently small A the spectrum is continuous. Given that kjWnj is a smooth 
function of aj for a £x n (as we numerically checked), we can apply the Riemann- 
Lebesgue lemma to extract the large time behaviour of (7.21). We hnd that it relaxes 
to 

Vnioo) ^Vn + ig-g) ^[A~%jVj . (7.25) 

j 

The power-law corrections to this result can be obtained by stationary phase 
approximat ions. 

On the other hand, when ±A exceeds some critical value ±A^, R’s spectrum 
develops an isolated eigenvalue and the corresponding eigenvector is bounded, whereas 
the continuous part of the spectrum is still described by (7.24). The main consequence is 
that the oscillations associated with the isolated eigenvalue do not cancel by dephasing 
mechanisms and local degrees of freedom keep oscillating at arbitrarily large times. 

7.1.2. The bound state. The isolated eigenvalue og can be easily worked out in 
momentum space 

AfnWnfi = aoWifl => [og - 16(1 + 5'g “ 2go cos k)]w{k; 0) = Ah(/c)wg,o , (7.26) 
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with w{k;j) = wqj + 2 cos(nfc)t(;„j and v{k) = no + 2 cos(nfc)n„. Isolated 
eigenvalues are such that to = y/a^/‘^ is outside of the continuous band, given by the 
image of \Jl + — 2go cos k. We are therefore allowed to write 

w{k]0) X v{k) 


Wo,o 16 — 1 — 5^0 + 2(70 cos k ’ 

which, integrated over fc, gives the condition 
^ _ A dk v{k) 

~ 16 


(7.27) 


(7.28) 


l_^ 27r — 1 — 5^0 + 2(7o cos k 

The critical values at the boundaries of the relaxation region correspond to the 
extrema of the continuous band where to = 1 + sgn(A)| 5 'o|- Using the explicit form of Vn 
we then find 

sin^ k 



^_ 

l-n 27r (sgn(A)|5(o| + 5-0 cosk)^/l + g^ -2go cosk^ 
and for A > A+ or A < A 7 equation (7.28) can be rewritten as 

r dk sin^ k 1 

J-TT 27r (o;^ — 1 — ( 7 o + 25 fo cos k) \/l + g^ — 2go cos k 2A 
Taking the derivative of this expression with respect to to we obtain 




0 < a; < |1 - l^foll V a; > 1 + IfifQ 


(7.29) 


(7.30) 


(7.31) 


this fact, together with the observation that A > 0 if a; > 1 + l^fol and A < 0 if 
0 < bj < |1 — l^foll; implies that A is an increasing function of uo (in the allowed 
dominion), therefore it is injective. This means that there can not he more than one 
isolated eigenvalue for a fixed A. 

The right eigenvector corresponding to the isolated eigenvalue is given by 
f" dk sin^fccos(? 7 ,fc) 


Wnfl oc 


(7.32) 


l-n 27r (tjS _ X _ ^2 2go cos k) ^/l + g‘^- 2go cosk ’ 
which decays exponentially with n, conhrming that it is a bound state. Analogously, the 
left eigenvector reads as (the factor in front of the integral is due to A’s asymmetry) 

dk cosink) 


w, 


0,n 


OC (2 ^no) 


27r — 1 — (yfg + 2gQ cos k 


n^—nO 


(2-M(-irc 

V(1 + -a;2)2 -4 c/2 


(A > 0), 


where 


9 = arccosh 


V 2(70 


(7.33) 


(7.34) 


Using (7.21) and (7.23) we can compute the entire contribution of the isolated eigenvalue 
to the solution of (7.19): 

V ■ Wq Wo 1 ~ cos(4a;f) 


(g-g)- 


wq ■ yjQ 


16a;2 


(7.35) 
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Figure 11. The time-dependent magnetic field h{t) (minus its long time average h) 
of the mean-field Hamiltonian (7.2) after a small quench from the ground state of 
the TFIC (7.13), with go = 1.5, and Hamiltonian H{g,X) (7.1), with g = 7.6513 and 
A = 3.5. The agreement with a function of the form h{t) ^ h ccos{£t + (p) becomes 
excellent at large times. The parameters of the fit are given by c « 0.000958(0.000989], 
E « 10.0507(10.0353], p = 0.0019(0], where in square brackets we reported the 
corresponding values based on the prediction (7.35). Despite the parameters do not 
differ much from the (asymptotic) prediction, at large times (inset) the corrections to 
the frequency have conspicuous effects. 


We stress that, within the linear approximation, the oscillation freqnency is 
independent of g. Fignre 11 shows that the most important correction to (7.35) lies 
precisely in the frequency, essentially because it is multiplied by the time, which has to 
be large for the subleading (time dependent) contributions to be negligible. However, in 
not-too-large time windows, the numerical data are in excellent agreement with (7.35) 
(the expression must be modihed including a corrective phase shift if {g — g)Jt is not 
small). 

We point out that there is a third relevant point A* < A7(< 0) at which a; = 0: 
dk sin^ k 


A* 


TT {l + gl- 2^0 cos/c)3/2 


-1 


(7.36) 


For A < A*, the isolated eigenvalue of A becomes negative [i.e. oj becomes purely 
imaginary). This would result in an exponential growth of (7.21), which after some 
hnite time would no longer be consistent with the linearisation procedure. However, 
we point out that, for small quenches in which the energy is close to the ground state 
one, A is always larger than A*. Indeed, if we assume A < A*(( 70 ), h turns out that 
the ground state of (7.13) is not equivalent to the ground state of H{g{gQ, A), A), with 
g{go, A) given by (7.17). The latter state is instead associated with the ground state of 
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the TFIC Hamiltonian (7.13) with magnetic held Qq ^ go such that g{g'o, A) = g{go^ A) 
and A*(( 7 o) < A. In the quench dephasing diagram of Figure 2 we can indeed easily 
identify A^ (the ‘critical’ curves for A positive and negative, respectively) but there is 
no trace of A*. 


Interpretation. The bound state of the matrix A may be put in relation with the 
existence of localised excitations in (7.1). We emphasise that this is not an ab initio 
calculation but rather a physical picture that explains the observations. 

Since the time evolution of the expectation value of any local observable in I'lpgg) 
is stationary, we can assume that, to all intents and purposes, \'ipgo) is an eigenstate of 
H{g{go, A), A). We now consider the limit of small quench and assume |A| > |Ac|. From 
(7.21) it follows that the projection on the bound state of y”{t)4:iuy'^{t) is proportional 
to an oscillating phase 

'^kniVnit) + • ( 7 - 37 ) 


The left hand side can be written as the expectation value ('0go(^)l-®gol’^9o(^)) ^ 

noninteracting operator with symbol 

{h{t) — cos{k))a^ — sin ka^ + iua^ 


sin k- 


U‘ 


-I- gl + 2go cos k 


(7.38) 


where go approximately satisfyes (7.16). In the no-quench limit the mean-held parameter 
is constant h{t) —)■ go (and go), so such operator can be written as 


i,n 


a| 



(7.39) 


where 




dfc Aqo — cos(k))a^ — sinfccr^ -|- iua^ 

—sm ^ ^- 

zn — 1 — g^ -\- 2go cos k 


(7.40) 


and a“ are the Majorana fermions (4.3). In (7.39) we left out the normalisation. 
Importantly, ~ with 9 dehned in (7.34); we have therefore found a quasi¬ 

local operator whose expectation value approximately oscillates in time as in (7.37). As 
a consequence, acts like an excitation over the initial state. Indeed we have 


{i’am'Ji’iM = E ~(7-41) 

£,n 


and if the state AA ^ sufficiently general representation (he. the overlaps with 
the eigenstates of H{g, A) are generally nonzero), (7.41) tells us that connects only 
states with energy difference equal to dcu. 

In conclusion, the bound state of (7.22) seems to be a manifestation of a localised 
excitation of H{g, A). 
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7.1.3. Remark on quenches from the ordered phase. In the previous section we ignored 
a subtlety that in principle could have invalidated part of the discussion (and part 
of the diagram in Figure 2). The mapping to a mean-held Hamiltonian relies on 
cluster decomposition properties but the ground state of the TFIC Hamiltonian in the 
ferromagnetic phase is the superposition of two Slater determinants [21], that separately 
do not possess cluster decomposition properties. Nevertheless, the discussion (and in 
turn Figure 2) remains correct also in this problematic case, at least for the operators 
commuting with This can be seen as follows: 

(i) The mean-held mapping is exact for the true ground state, which breaks the spin- 
hip symmetry realised by ‘^1- Thus (7.3) is valid. 

(ii) Using that is a quadratic operator in the Jordan-Wigner fermions, in (7.3) the 
ground state can be replaced by one of the two Slater determinants, therefore 
h{t) is still solution of (7.7). Analogously, the expectation value of any operator 
commuting with cr| can be found replacing the initial state with one of the two 
Slater determinants and then using Wick theorem. 

(iii) For any given t, the expectation value of operators 0°(£) that anticommute 
with Yl^cr^, like the order parameter, can be obtained from the large r limit of 
(To|t/^F(f)0°(£)0'’(£ + r)UMF(t)|To), using cluster decomposition properties. 

However, using similar general arguments we are not able to exclude that the expectation 
values of the odd operators might keep oscillating also when all the even operators 
relax. There is indeed a subtle problem of limits that comes from the trick of Point 
(iii). Nevertheless, in the cases considered we have never encountered this situation, 
suggesting that for the model under examination such complications do not arise. Here 
we provide a heuristic argument. When the limit of inhnite time for h(t) exists, at 
sufficiently long times the dynamics is essentially the same as for a quantum quench 
in the TFIC from a certain state with cluster decomposition properties. In the latter 
situation we can apply (a direct generalisation of) the results of [19], which showed that 
the expectation values of odd operators decay to zero. This suggests that the diagram 
of Figure 2 could be valid also taking into account the odd operators. 

7.2. Relaxation properties 

We now focus on the regions (of the parameter space) in which the limit limi^oo h{t) = 
hoo exists and investigate more closely the relaxation properties. Even beyond the linear 
approximation (7.21), we can still guess an asymptotic form for the solution of (7.7) 

2 /fc Cfc + -|- h.c.), (7.42) 

where Sk 2^/1 + — 2hoo cos k. Equation (7.42) is compatible with the relaxation of 

h(t) with corrections 0(t~^^‘^~^) (stationary phase approximation), where j is an integer 
that depends on the behaviour of Ak around the extremal points of the dispersion 
relation {k = 0 V tt). Our numerical analysis for several quenches from TFIC initial 
states is compatible with j = 1. This is not surprising since the same exponent governs 
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Figure 12. The time-dependent magnetic field h{t) of the mean-field Hamiltonian 
(7.2) after a quench from the ground state of the TFIC (7.13), with go = 1.5, and 
Hamiltonian H{g,X) (7.1), with g = 0.5 and A = 0.5. The dashed red line is the 
asymptotic prediction (7.43) with the coefficients estimated by fitting the data at 
very large times t > 1800. The inset shows the goodness of the prediction in some 
intermediate time window. The horizontal grey line corresponds to the value of hoc 
extracted from the fit. 


the late time behaviour of (a^) after quenches in the TFIC [21]. Thus, we conjecture 
the large time expansion 


h{t) 


~ hoo + 


Ao cos(4(l - hoo)t + (fo) + A^ cos(4(l -h h^)t + ip^) 


(7.43) 


and similar behaviours for yn{t). Remarkably, the (leading) oscillatory frequency is only 
determined by h^. In Figure 12, (7.43) is compared against numerical data for a quench 
that leads to relaxation. 

Since for asymptotically large times the time evolution is equivalent to the one 
generated by the TFIC Hamiltonian 


i 


(7.44) 


at late times local observables are described by a generalised Gibbs ensemble constructed 
with the local conservation laws of 77/. However, the Lagrange multipliers can not be 
simply hxed by computing the corresponding integrals of motion at the initial time, as 
they are in fact conserved only at asymptotically large times. This is an example of 
a stationary state written in terms of operators commuting with one another but not 
with the Hamiltonian. In fact, this is not the hrst time that such an unusual description 
emerges: in [61] the pre-thermalisation plateau was described by a GGE constructed 
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with operators in involution that are however not conserved at the perturbative order 
that was worked out. 

As it will be clarified in the next section, such a stationary state coincides with 
the pre-thermalisation plateau of [60]. In the thermodynamic limit this is just the 
stationary state that emerges at infinite time after the quench. Indeed, in the regions 
in which there is relaxation, we have not found any indication of pre-relaxation/pre- 
thermalisation behaviour. Furthermore, ‘relaxation’ is not synonym of ‘thermalisation’. 
Indeed it is not difficult to show that at late times the system still retains infinite 
information about the initial state. To this aim, as initial state we choose the ground 
state of the local Hamiltonian 

n 

H[{a}\ = Y.a, Hj , Oj G M, (7-45) 

3 

where n is finite and Hj are the most local reflection symmetric conservation laws of a 
TFIC model with Hamiltonian Hq. Such state can be easily constructed [82] and, using 
the notations of Appendix E, is completely characterised by the function 


m(/c) = —sgn(^^ Oj cos(j/c)), (7.46) 

3 

which is equal to 1 if and only if the excitation a\ of Hq is present in the state. It 
is important to note that different characteristic functions m{k) correspond to locally 
inequivalent states. Thermal-like behaviour would imply that the only information 
about m{k) that is retained at late times is the corresponding energy and magnetisation. 
We now consider the special cases in which the initial magnetisation is such that Hq 
is the mean-field Hamiltonian at the initial time. In this way, by Corollary 5.5, the 
expectation value of local observables is independent of time and the late time stationary 
state is equivalent to the initial state. The only scenario compatible with thermal-like 
behaviour is that each distinct function m{k) of the form (7.46) corresponds to a distinct 
pair where hm is the parameter of the late time mean-field Hamiltonian and 

Em the energy. The self-consistent conditions behind the no-quench limit are worked 
out in Appendix E and are given by 


hm — S' + 2A 


^m. — 




9 


4A 


dfc , . hm —cos k 
27r ^yl + hl^- 2hm cos/c 
dfc hm cos k — 1 

..27r 


(7.47) 


(7.48) 


2hm cos k 

It is easy to see that the same pair {h-m, Cm} is associated with infinitely many functions 
m{k) of the form (7.46) (it is enough to choose n = 3 in (7.46) to find (infinite) 
examples). Thus, the non-equilibrium time evolution under (7.1) does not generally 
result in thermalisation. 

We conclude the analysis of (7.1) considering the expectation value {Qn)t of the 
local conservation laws of 77/ (7.44), which can be written as follows [19] 


Qn 

~T 


^ cos{nk)£k (aia, 

ZTT \ 



(7.49) 
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where Sk is the dispersion relation of if/ and Uk are noninteracting fermions that 
diagonalise Hf. Using free-fermion techniques we obtain 

g0n. 

In particular, the expectation value of Hj = Qq can be written as follows 


/Qn \ ^oo ~ 


(7.60) 


Pf\ 

i-r) 

L t 


{h{t) 


h U 

"■oo j 


[hoc -gf 


4A 


(7.51) 


from which it is clear that the relaxation exponent of (if/)^ is twice the exponent of h 
(which in the cases that we investigated is 3/2, cf. (7.43)). This result can be easily 
generalised to any local conservation law of if/. By taking the time derivative of (7.50), 
the last equation of (7.10) implies 

' (7.52) 


/Qn\ _ h 

'Xu ^ 2 


'OO / 


Since both h — and y'^ are expected to decay as t (with oscillatory factors like 
in (7.43)), we immediately obtain 


(%> - {%) ~0(« 

L/ oo ij t 


-3\ 


(7.53) 


7.2.1. Comparison with [60]. It is not a coincidence that the same relaxation exponents 
were found in [60] in a perturbative framework. In order to understand the relation 
between the two models we must come back to the modified version of Hamiltonian 
considered in [60], i.e. (2.8). The mean-field Hamiltonian for that precise model reads 

HufH) = - ^ - g ^ cr| + 2A (cr| - a] - ^ cr|) , (7.54) 

i i. i i 

where the time average ■ • ■ is taken with respect to the Hamiltonian with A = 0. It is 
convenient to introduce the auxiliary Hamiltonian 

Hupit) = — — g 2A {a^ — jyjp ~ ’ (7.55) 

i i i e 

where the time average ■ ■ ■ is now taken with respect to Hf = lim^^oo HMF{t), limit 
which is assumed to exist. It is now simple to prove that Hf = H{g, 0), where H{g, A) 
is given in eq. (2.8). To this aim let us consider the expectation value of evolving 
via ifMF(t); it fulfils 

t—KXi ’ t—^co T—>-oo I Jq 

= lim lim — / ds {'^o\Uljp{ts)a]UMF(t + s)\^o) 

t—>-oo T^oo 1 Jq 

= lini = (V) ]y[p , (7.56) 

t^oo ’ ’ 

where UMFif) is the time evolution operator constructed with HMF(t) and in the second 
step we replaced e~^^f^UMF(t) with UMF(t + s), as it is legitimate at late times. From 
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this it follows Hf = H{g,0) and, in turn, the equivalence between (7.54) and (7.55) 
77 mf(^) = 7/mf(^)- Importantly, this means that (7.54) and (7.13) have the same inhnite 
time limit if we set hoo = g- 

As a matter of fact, the equivalence between (7.54) and (7.13) is not restricted to 
inhnite times. The mean-held time evolution operator for (7.54) can indeed be written 
as follows 

UuFit) = Hf + 4X S%s)^ . (7.57) 

Here ^ ct|, 6a^ = and 

S^(s) = _ (7.58) 

If the magnetisation relaxes faster than l/f^’*'", with a > 0, the operator at the exponent 
of the hrst term of (7.57) is a bounded function of the time, so that exponential can 
be safely expanded. In fact the entire term can be neglected (it gives corrections o(A)). 
The same holds true in 5*^(5), indeed the hniteness of sup^ A ds' {Sa^) j^p guarantees 
that the series expansion of the exponentials in S^{s) can be truncated for any s with 
an error that goes to zero as A —?• 0. By considering the hrst terms of the expansion 
one immediately realizes that the correction is o(A) and approaches zero for large s as 
1/s". Putting everything together, replacing S'^(s) by in (7.57) introduces an error 
o(A), independently of the time. Therefore we obtain 

f7MF(^) ~ Texp^-i y ds 77/4A (5(7^)^ j^p -h o(A). (7.59) 

We can also replace {5 o'^)smf with (cr^) 5 MF ~ (‘^^)oo mf- diherence between the two 
terms is o(A°) and approaches zero at large times at least as In this way we have 

reduced the Hamiltonian (7.54) to (7.1), provided that the condition (see hrst equation 
of (7.10)) g = g + 4Am^ is satished. With this choice we recover the perturbative results 
of [60], e.g. the relaxation exponents (7.53). A posteriori we note that the large time 
behaviour of (ct|) under (7.1) {cf. (7.43)) is sufficiently fast (a = 1/2) to justify the 
approximation of (7.57) by (7.59). We also checked that the mean-held solution of (7.1) 
is perfectly compatible with the results shown in Figure 1 of [60] (see Figure 13). 

Finally, we point out that [60] introduced the term some conditions in 

the long-time limit, where (7.1) and the Hamiltonian of [60] turn out to be equivalent. 
Equation (7.1) is therefore a sensible replacement for the Hamiltonian of [60]. 

1.3. Generalisations 

The construction of (low-entangled) stationary states that we proposed in the previous 
section ((7.45) and below) can be applied also to other Hamiltonians of the form (4.9) 
if the corresponding mean-held Hamiltonian is integrable. In those cases we can rule 
out thermalisation if the self-consistent problem satished by the stationary solutions at 
hxed energy and mean-held parameters has more than one solution. 
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Figure 13. The time evolution of the number of quasiparticles Uqp = that 

diagonalise the late-time mean-field Hamiltonian for a quench with go = 8, X = 0.05 
and g = 3.597274.... The parameters are chosen to reproduce the first figure of [60] 
(dashed orange line). In particular, g is such that the mean-field parameter h(t) in the 
limit of infinite time approaches the value g = 3.5, considered in [60] (see the main 
text). The timescale and A differ from [60] because of two small typos (the dispersion 
relation was unintentionally halved and the right hand side of Equation (3) of [60] 
should have been multiplied by 4). The (tiny) discrepancy is compatible with higher 
order corrections in A. 


When the mean-field Hamiltonian is noninteracting, following the lines of the proof 
sketched for the nonlocal generalisation of the Ising model in Section 7.2, one can 
generally show that the solution is not unique. 

It is also reasonable to expect that, also in the presence of interactions, the finite 
number of constraints given by the energy conservation and the late time values of 
the mean-field coupling constants could not reduce the parameter space of the initial 
Hamiltonian (7.45) to a single point. We indeed believe that thermalisation is unlikely 
to emerge if the mean-field Hamiltonian describes an integrable model at any time. 

Nevertheless, the interacting case exhibits counterintuitive behaviours, for example 
in the energy level-spacing statistics. Generally integrable models exhibit Poisson 
statistics, whereas generic models follow a Wigner distribution [83, 84]. There are 
many exceptions to this rule [85] , however the nearest neighbour spacing distribution is 
probably the most reliable numerical check of integrability. 

In Figure 14 the level spacing distribution is shown for various chain sizes for the 
Hamiltonian (7.1) with A = 0.5 and g = 0.5 in the reflection symmetric sector of the zero 
momentum subspace with spin-flip parity equal to 1. The numerical data suggest 

that in the thermodynamic limit the curves collapse to an exponential distribution 
(Poisson statistics). This is consistent with our observation that at arbitrarily large 
times after a quantum quench the system keeps retaining infinite information about the 
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Figure 14. The nearest neighbour spacing distribution P{s) of the Hamiltonian 
(7.1) with g = 0.5 and A = 0.5 in the reflection symmetric sector of the zero 
momentum subspace with spin-flip parity ct| equal to 1 for various chain sizes 
(L = 10 19). As the size is increased the colours vary from green to brown and 

the lines become thicker. Dashed lines correspond to odd sizes. The dotted black 
line is the exponential distribution (Poisson statistics). In the inset the distance 

d{P,Pstat) = \Jfi^(P(s) - Pstat{s)y from Poisson (Pstatis) = e~^) and Wigner 
(Pstatis) = fse“T'*^). 


initial state. 

A completely different scenario appears for the Hamiltonian 


H = 




1 + 7 




+ 


1-7 


^I 


A 

^ 4L 






(7.60) 


^ i 

The corresponding mean-field Hamiltonian for a given one-site shift invariant initial 
state I To) is given by 


+ V^77+i) +1 ■ (7.61) 

i e 

At hxed time, this describes an XYZ model, which is known to be integrable for 
any choice of the coupling constants. Therefore, assuming relaxation, the stationary 
properties of local observables should be described by a GGE constructed with the 
conservation laws of an XYZ model. We considered 7 = 0.25 and A = 0.5. We found 
large hnite size corrections in the level-spacing statistics (in the same sector as before) 
and, in particular, a remarkable even-odd parity effect (see Figure 15). However, it 
seems that increasing L the curves collapse to a Wigner distribution, as it commonly 
happens in non-integrable models. 
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Figure 15. The nearest neighbour spacing distribution P{s) of the Hamiltonian (7.60) 
with 7 = 0.25 and A = 0.5 in the reflection symmetric sector of the zero momentum 
subspace with spin-flip parity equal to 1 for various chain sizes (L = 10 19). 

As the size is increased the colours vary from green to brown and the lines become 
thicker. Dashed lines correspond to odd sizes. The dotted black line is the Wigner 



distribution. In the inset the distance d{P,Pstat) = y iP{s) — Pstatis))'^ from 



Poisson (Pstat(s) = e“®) and Wigner (Pstat(s) = f se“T®^). The distribution for odd 
chains converges rather quickly to Wigner. 


This might appear in contradiction with our conjecture that thermalisation should 
not be expected when the mean-held Hamiltonian is interacting and integrable. In fact, 
we have not taken into account that at late times the mean-held parameters are hxed. 
In the previous section we ruled out thermalisation by constructing an inhnite family 
of stationary states with the same energy and the same mean-held parameters. If this 
is possible, then we should hud a signature of the huge degeneracy in the level-spacing 
statistics by restricting the space to the excited states that lie in some shell with mean- 
held parameters almost hxed. 

In the restricted space our preliminary analysis is indeed compatible with Poisson 
statistics also for the Hamiltonian (7.60). However, our data turn out to be compatible 
with Poisson statistics even if the mean-held Hamiltonian does not describe an integrable 
model. This is in contrast to our expectations that in generic systems there should not 
be more than a few parameters that characterise the stationary state, namely the energy 
and, at worst, the mean-held coupling constants at inhnite time after the quench. 

Our interpretation of these contradictory results is that we did not investigate 
sufficiently large chains, so our analysis of the energy-level statistics in the restricted 
space is not sufficiently indicative. We are conhdent that a more accurate analysis will 
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show a different behaviour in the non-integrable case. 

Finally, we point out that the situation is trickier when there are isolated points in 
the parameter space of the mean-field coupling constants that correspond to integrable 
models. For example, it is not clear to us whether or not we should expect thermalisation 
when at asymptotically large times the coupling constants of the mean-field Hamiltonian 
match the integrability points. 

7 . 4 . Summary 

We showed that the time evolution under (7.1) has a quite rich phenomenology, including 
both cases of relaxation and cases of persistent oscillatory behaviour. In the limit of 
small quench the latter has been interpreted as the effect of localised excitations that 
appear (or become relevant) when the Hamiltonian parameters cross some “critical” 
line. 

In addition, we confirmed the perturbative results of [60] in a non-perturbative 
setup. Our analysis excludes that in the thermodynamic limit the late time behaviour 
of local observables could be described by a thermal-like ensemble. 

More generally, we provided some argument that suggests that thermalisation is 
unlikely to emerge if the mean-field Hamiltonian describes an integrable model at any 
time after the quench. 

We also proposed a numerical check of thermalisation based on the analysis of the 
energy-level statistics on some restricted space, however our preliminary analysis on 
small chains (L < 20) was not sufficient to discriminate between the cases in which we 
expect thermal-like behaviour and the cases in which instead also at late times infinite 
information about the initial state is retained. 

8. Conclusions 

Pre-relaxation is a dynamical phenomenon that arises when small perturbations break 
symmetries that affect the late time behaviour of local observables. When the 
perturbation breaks (abelian) integrability, this is usually called pre-thermalisation, 
which is generally thought as a two-step process in which local observables experience 
virtual relaxation before approaching thermal-like expectation values. However the 
relaxation process can also be more complicated, following many steps of quasi¬ 
stationary behaviour. This happens in particular when the model is close to a non- 
abelian integrable point. In order to extract the pre-relaxation behaviour one must 
therefore identify the correct time scale of the phenomenon. 

We have considered the problem of pre-relaxation after quantum quenches in weakly 
interacting models, starting from initial states with cluster decomposition properties. 
We focussed on the particular situation in which the unperturbed (one-site shift 
invariant) Hamiltonian has a non-abelian set of local conservation laws that break 
one-site shift invariance. In particular we considered interacting perturbations to the 
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XY spin-^ chain and investigated both integrable extensions, like the Heisenberg XYZ 
model, and the effects of perturbations that break integrability. 

We identified the inverse perturbation strength as a relevant time scale of pre¬ 
relaxation and studied the dynamics of local observables at times proportional to it. 

Despite the model being interacting, the noninteracting structure, remnant of 
the unperturbed Hamiltonian and manifested in the Wick theorem, survives the pre¬ 
relaxation limit. However interactions do affect the dynamics by introducing a nontrivial 
time dependence in the effective noninteracting Hamiltonian that generates the time 
evolution. The most striking effect is probably that, even if local degrees of freedom 
approach stationary values, these can not be generally predicted without following the 
entire dynamics. 

We have shown how to recast the non-equilibrium problem into a system of 
nonlinear differential equations involving expectation values of quasi-local operators. 
The system of equations has qualitatively distinct solutions, which vary from trivial 
stationarity to persistent oscillatory behaviour over the entire time window considered. 
We have not found any relevant difference between integrable and non-integrable 
perturbations, suggesting that the scenario of thermalisation in generic models arises at 
much larger times. 

For the very nature of the local conservation laws of the XY model, in order to 
have a nontrivial time evolution the initial state must break one-site shift invariance. 
For a particular initial state of that kind we considered a limit in which the equations 
can be linearised and exhibited the analytic solution, in which one-site shift invariance 
is eventually restored. The regime worked out analytically shows quite clearly the 
importance of cluster decomposition in the non-equilibrium problem. While, as 
mentioned above, the pre-relaxation limit is trivial for one-site shift invariant states, 
a shift symmetrisation of the two-site shift invariant initial state has a nontrivial time 
evolution. This is because cluster decomposition has been lost with the symmetrisation. 
It is important to take into account such aspect when analytic predictions of late time 
stationary behaviour are compared with numerical data at times in which one-site shift 
invariance is not yet restored. 

The crossover between oscillatory behaviour and relaxation is quite interesting per 
se. This has been the main motivation for the analysis of a simplified model that 
shares most of the formal aspects with the effective description of pre-relaxation in the 
perturbed XY model, but that, in fact, has not been derived from a pre-relaxation limit. 
We considered a transverse-field Ising chain with an additional nonlocal interaction 
proportional to the magnetisation squared per unit length. This model was already 
studied in [60] in the framework of a perturbation theory. We used some general 
properties, proven for Hamiltonians of that form, to obtain nonperturbative results and 
showed that in the thermodynamic limit subsystems retain infinite information about 
the initial state, whatever large the time is. This is not in disagreement with [60], where 
thermalisation was conjectured for time averages in finite systems: for this model the 
diagonal ensemble could not be locally equivalent to the stationary state that emerges in 
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the thermodynamic limit when the quench parameters are compatible with relaxation. 


We showed that the late time behaviour in the thermodynamic limit (which 
corresponds to the ‘pre-thermah behaviour of [60]) can not always be described by 
a stationary state. In the parameter space there are indeed ‘critical lines’ that separate 
relaxation from persistent oscillatory behaviour. We dehned a limit of small quench and, 
in that limit, exhibited the analytic expressions for such critical lines. The appearance of 
oscillatory behaviour has been interpreted as a consequence of the emergence of localised 
excitations. 

We also discussed the generalisations to other Hamiltonians in which some terms 
have the form of interactions with macroscopic observables, like the magnetisation 
squared per unit length of the model above. In particular, we ruled out thermal-Iike 
behaviour in a large class of models of that kind. 

Finally, we would like to stress that our description of the pre-relaxation limit 
is based on a few hypotheses. In particular, we neglected some “anomalous terms”, 
proving only the self-consistency of the conjecture. Some preliminary checks against 
iTEBD simulations are conhrming the validity of the assumptions [89]; however, we 
leave a more rigorous analysis of the regimes of validity of our approximations to future 
research. 
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APPENDICES 
A. Free-fermion relations 

We briefly summarise some useful relations valid in noninteracting models. Additional 
details can be found in [65] and [86]. 

Let us consider a local one-site shift invariant spin-chain Hamiltonian H that is 
mapped to noninteracting fermions by the Jordan-Wigner transformation ({a", a^} = 





(A.l) 


j<e 
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Ljn 


^ ^nm—n+1 
nm—n+1 




^nl-n+l ®n£-n+l ••• ^nl ^]^'rn 


£,m 


where L is the chain length and "H is a block-circulant^ matrix 

[H^%m =jYl (k) = 1 , 


where 




^nm 

\ <m ) 


, (A.2) 


(A.3) 


(A.4) 


The index n is a divisor of L (but in the thermodynamic limit any positive integer is 
allowed). We call the 2?7,-by-2n matrix 'H^'^l{k) the n-site representation of the symbol, 
which completely characterises the block-circulant matrix. In the following we will refer 
to also as the symbol of H. 

For given n, the following properties hold: 

(i) Any function of block-circulant matrices is a block circulant matrix, with symbol 
equal to the function of the respective symbols 


f{A,B,...)^f{A{k),B{k),...). 


(A.5) 


(ii) Let A and B as in (A.2). Their commutator [A, B] has the form (A.2), with symbol 
equal to the commutators of the symbols 


[A,B]^[A{k),B(k)] 


(A.6) 


(hi) The time evolution in the Heisenberg picture under (A.2) of a noninteracting 
operator A of the same form is noninteracting, with symbol 

^iHtAe-iHt (A.7) 


(iv) (Wick theorem) The expectation value of any operator in a thermal state of (A.2) 
(and in any Slater determinant state) can be expressed in terms of the correlation 
matrix 


p(n) 

^ Im 


/ 

^nm— 


Siml - ( 


nm—n+1 ^ 

V 
1^ 

^nm—n+1 


/ 


'^n^-n+l 


<l )) • (A.8) 


ji A circulant matrix M is a Toeplitz matrix = Mn-m ) in which any row is a right cyclic shift 

of the row above. 
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(v) For a thermal state with inverse temperature (3 the correlation matrix is given by 

= — tanh(^'H). (A.9) 

Therefore, the “thermal” ground state (which, in the presence of degeneracies, is 
equivalent to the incoherent superposition of the states) has correlation matrix 

Too = -sgn('H). (A. 10) 

(vi) If Pi is block-circulant, F is block-circulant as well, so it is completely characterised 
by its symbol F(fc). In particular, the symbol of the thermal ground state reads 

Too{k) =-sg\i{PL{k)). (A.ll) 

(vii) The correlation matrix of a Slater-determinant state that time evolves under the 
noninteracting Hamiltonian (A.2) has the symbol 

F(fc; t) = 0)6*^^^^*. (A.12) 

viii) In the thermodynamic limit {L —)■ cx)), the expectation value of an operator A 
of the form (A.2) in a Slater-determinant state with correlation matrix F can be 
written as follows 

1 1 dh 


B. Time averages of interacting operators 


In this appendix we show the validity of Property 4.1. The Property can be more easily 
proven for Jordan-Wigner fermions 

(BA) 

the relation for the Majorana fermions aj, aj will then follow by linearity. 

In order to proceed it is convenient to introduce the following notation 


c“(t) 


c](t) a = + 
Cj{t) a = — . 


(B.2) 


The relation between c\c and the Bogoliubov fermions b^{k),b{k) that diagonalise the 
unperturbed (noninteracting) Hamiltonian Hxy (4.1) can be written as 


c“(t) 


1 


k 0=± 


(B.3) 


Here U{k) is the 2x2 matrix dehning the Bogoliubov transformation, Sk is the dispersion 
relation in the 1-site shift invariant representation 


J 


\J cos^ fc -|- 7^ sin^ k 


and we set 

r W{k) 13 = + 

I K-k) (3 = -. 



(B.4) 

(B.5) 
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The relation (4.6) is then eqnivalent to 

j^Z^'-j+nl'-j+n2'-j+n3'-j+ni “ {n} ^ '^{n} ’ 

j 

where is a factorised term 

{n| 

-r{“} - V 

-^{n} - 

s=0 

<A = \ '^{-^Y^^"+nA+n, (B- 8 ) 

s ^ 

and is the remaining contribntion. Using (B.3) we can explicitly carry ont the 
time average and the snm over ^ in (B. 8 ). We obtain 

S ^ 

where a = —a and we defined 


Ql 0:2 p03 04 

^71x712 ^7lZ^7l4, 


^Oi _03 I 

'^ni'^n3 '^712 ^^^4 ~'~'^n2'^n3 '^ni'^n4 


., 02^03 ^01^04 


(B.7) 


kg = k + ns . 

Analogously, (B. 6 ) reads as 
1 


E- 


„«! „02 „«3 p«4 

^ i+«i i+?i 2 7+^3 i+”-4 


(B.IO) 


= E n bft(fa)b&(fe)bfe(A^ 3 )b^,(fa) 

{ki} j=l 

X(5oifci+a2A:2+03fc3+04fc4'^/3iefcj+/32£fc2+/33£fe3+^4efc4 ' (B.ll) 

In order to compute the sums over the momenta it is necessary to solve the constraints 
given by the delta functions, i.e. 


C(iki T 02^2 T T 04 / 1:4 — 0 


YlSki + /^2^fc2 + + /^UA;4 — 0 . (B-12) 

Some of the solutions to these equations can be found by requiring the terms of (B.12) 
to cancel in pairs. This would give 

'5aifcl+a2A:2+«3fc3+«4fc4*^/3l£fcj+;32£fc2+/33£i;3+/34£fc4 ^ ^ ‘^1 ~k ^2 A Ag , (B.13) 

£ = 0,1 


with 


^ (Si,^2^133,y4,^^lki,aik2,s^a3k3,a4,k4^s (®A4) 

A 2 *^/3l„S3*^/32,A*^"l^l'“3fc3,s*^<52fc2,04^4,3 (B.15) 

A 3 ^/3i„a4^/32„33^“lfcli«4fc4,3^(52fc2,03^3,8 ' (B.16) 

For a generic dispersion relation it is reasonable to expect these solutions to be the only 
possible. For the specific dispersion considered, (B.4), equations (B.12) admit other 
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solutions in the thermodynamic limit. We call these anomalous solutions because they 
strongly depend on the precise form of the dispersion relation. We now show that the 
term in (B. 6 ) is exactly the contribution arising from these solutions, i.e. 

, 4 

Pn) = 15 E E n hf,{k,)h^Ah)hi>,(h)hf,(k,) 

ki,k 2 ,k:^,k 4 s=0,l j=l 

x{Al + A^ + A^}. (B.17) 

We stress that the operator will be nonzero only in the thermodynamic limit. 

We consider for example the term containing A 2 . We have 


A:i,/c2,fc3,/c4 s=0,l j=l 

p,q s=0,l 

X (p)b^2 (g)b^^ {aia3ps)h^^ {a2a4qs) 


E 

s=0,l 


Oil as 012^0:4 
'^ri.2'^n4 


+ 0(L- 


(B.18) 


where we used the commutation relations of the {hp{k)} in the last step. Although 
the terms and are in fact multiplied by L in the time average of (4.2), the 
possible corrections 0{L~^) in (B.18) (which would result in corrections 0{L^) in the 
effective Hamiltonian) are locally irrelevant, because their density approaches zero in 
the thermodynamic limit. 

We obtain analogous results for Ai and A 3 , that is to say (B. 6 ). 


Remark We point out that for other dispersion relations (still with the properties 
Sk = Sk+Tt and Sk 7 ^ Sk+ir/n for generic k and u > 1 ) the anomalous terms could be 
factorised as well. Generally in such situations the factors have a very simple time 
dependence, e.g. a single oscillation frequency. As a consequence, relaxation is ruled 
out. 


C. Towards a mean-field description 

In this appendix we prove the Lemmas of Section 5. 

Lemma C.l If O E S, the operator norm (^i.e. the maximal eigenvalue in absolute 
value) ofO/L is hounded. 

Proof The proof is straightforward. Let us expand OfL as in (4.9): 


N 


j=l 


(C.l) 


where Om have local densities, that is to say, they can be written as follows 

= (C- 2 ) 
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with local operators. We immediately find the chain of ineqnalities 


n ^ ^ ^ 


Ol- 


U) 


o 


U) 


N 


ml ’ 

The right hand side is clearly 0{L^) because N and Uj are finite by definition and 
are local. 

Lemma C.2 If 0,0 & £, then [O, O] E S as well. 

Proof Without loss of generality, we can restrict to two single terms of the expansions 
(4.9) of O and O. The commutator of the two terms reads as 

^ nO _ _ _ /nh) ^ 




i=i 


[-—-of ■ ■ ■ of, -—.of ■ ■ ■ of] = 

^ Tni — 1 ' J rij — l ^3 -J 


l^nj+nj-2 


Y, of ■ ■ ■ of_,6f ■ ■ ■ oHfiof, d«]o<!>i • • • of off ■ ■ ■ of. (c,4) 


k^p 


Since [Of ,Of] have local densities (the commutator of two local operators is nonzero 
only if there is a region on which they both act nontrivially; in addition, its range 
is smaller than the sum of the ranges of the two operators), the number of extensive 
operators exceeds by one the exponent of 1/L. Thus, [O, O] E £. 

Lemma C.3 (viz. Lemma 5.1) Let O E £ and |\h) a state with cluster decomposition 
properties. The expectation value ofO/L in |\I/) can be reduced to the expectation values 
of the local translation invariant operators it consists of: 


’’ L^oo • 


L —^oo 


(C.5) 


Proof Let us consider a term (5.1) of the expansion (4.9). Its expectation value (per 
unit of length) is given by 




(C.6) 




where hj^i- are local operators acting nontrivially only around £j and such that 

h,^Y ■ (C:-’’) 


By cluster decomposition we have 




(4-1- 






-/«,L),(C.8) 


iq-q/l>e»i (vjVjO 




where 


lim lim ff,L) = 0 . 

L^OO 


(C.9) 
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The difference between (C. 6 ) and the left hand side of (C. 8 ) can be bounded from above 
as follows 


/I<5 

Analogously 


-IT) 




< 

— max 

" V2y 

L {£} 


—^O.(C.IO) 




E 


($1(11,,J*) {SIAnAl*) 


. fn^ 

^ 1 

^ 1 „ 

— max 

" V2y 

'L ml 




i=i 


so that 




(Sl/fjl®) 


< |/«,L)I + 0 ( 1 /L) 


(C.12) 


Being ^ arbitrary, we can take the limit lim^_).oo lim£,_).oo, obtaining (C.5). 


Lemma C.4 If the state |To) has cluster decomposition properties and O & £, the 
mean-field Hamiltonian defined in section 5 satisfies the following identity: 

lim Tl£'(i)l»o> = lim {MUHOIW t1C^(*)I'1'o) , (C.13) 

L^oo 1j L^oo Jj 

where U was defined in (5.4). 


Proof Let us consider a generic term (5.1) of the expansion (4.9) of H 

H = -^Hi ■■■Hn. 
rn—1 ^ "■ 

The corresponding term (5.3) of the mean-field Hamiltonian (5.3) is given by 

= E n , 

£=i j^e 

By taking the commutators with O we find 




n £—1 


["Ti=5:n 

£=1 j=l 


Hj[H£,0] 

L L 


n 


n 

j=£+l 


L 


^ e=i ^ ^ 


(C.14) 


(C.15) 


(C.16) 

(C.17) 


Because |To) has cluster decomposition properties and the mean-field Hamiltonian is 
local at any time, the state U{t) |To) has cluster decomposition properties as well 
(the only difference with respect to |To) is that the function / of (C. 8 ) is now time 
dependent). Finally, by Lemma C.3, in the thermodynamic limit the expectation values 
of (C.16) and (C.17) in the state U{t) |To) are identical, that is to say (C.13). 
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Lemma C.5 Let |'&o) a translation invariant state with cluster decomposition 
properties and H,0 G £. The time derivatives of the expectation value of O/L in 
the state evolving with H^p{t) fulfil 

d"^ - O - _ O - 

-—\\m (>I>„|£/'(«)yf7(i)|>I'o>=*” lim (C.18) 

Ufi Lj - ^OO Jj L/ -^OO J_j 

n 

Proof We proceed by induction. First of all we see that for n = 0 the property is 
trivially satished; let then the property be true for n, we have 
d _ fO _ H - fO _ 

^ Jto {>I-„|et(i);^(7(«)l'I'o) = c- ^lim (<l-„|gt(i) |g,|g,...|g . -]...]|g(i)|<I-„) 

n 

= hm {^,\W{t)[HMF{t), [H, [H, ...[H, %.]]U{t)\^o) 

L^oo s^ L 

n 

= Hm ($o|g«(t) [g, |g, ...|g , g]...]]g(t)|'l.o) , (C.19) 

n 

In the second step we used Lemma C.2 and Lemma C.4. This concludes the proof. 


Lemma C .6 (viz. Lemma 5.2) Let |To) he a translation invariant state with cluster 
decomposition properties and H,0 G £. Let the expectation value of O in the state 
that time evolves with H^p{t) be an analytic function oft in the strip |Im[t]| < r, with 
r a nonzero constant. In the thermodynamic limit, the time evolution with H can be 
replaced by the time evolution with the mean-field Hamiltonian: 


lim 

L —^CXD 1^ 

= lim (To|t/^(t)^f/(t)|To) . 

L^OO Ij 

(C.20) 

Proof We dehne 



fit,s) = hm (vl/o|m(t)e*^^-e-‘^^t7(t)|v]/o) . 

L^oo ij 

(C.21) 

By Lemma C.5 we have 



Qn Qn 



^^n/(^’ 0 ) - 

o 

II 

(C.22) 

indeed 




(C.23) 


n 


By assumption, f{t, 0) (which corresponds to the time evolution with the mean-hied 
Hamiltonian) is analytic in the strip |Im[t]| < r, so the convergence radius of the Taylor 
expansion at t = 0 is larger than or equal to r. Thus we have 

Qn 


/(fO) = 


pTi Qn 

n\ dP 


/((.o) = E 

t=o ^^ 


r 


n\ dP 


t=o 


/(O,^) = /(0 ,r), 


T 


<r. (C.24) 


Let us call a time such that /(t,0) = /(0,t) for any 0 < t < t*. As before, the 
function f{t + r, 0) is analytic in the strip |Im[r]| < r, so we have 


fit + r, 0) = 0) = = /(O’ ^ ’ 


|r| < r . (C.25) 
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That is to say 

/(t,o) =/(o,t) yt<p ^ /(t,o) =/(o,t) yt<p + T. (c. 26 ) 

Since r is finite and (C.24) holds, we conclude 

/(t,0) = /(0,t) Wt, (C.27) 

which is exactly (C.20). 

Corollary C.7 (viz. Corollary 5.3) Lemma C.6 holds true in particular for local 
operators. 

Proof By translation invariance, the expectation value of any local operator O is equal 
to the expectation value per unit of length of the operator E £, obtained by shifting 
O along the chain and summing all the (L) terms. 

Corollary C.8 (viz. Corollary 5.4) Let |\I/o) a translation invariant state with cluster 
decomposition properties and H E £. In the thermodynamic limit, the time evolution 
of the reduced density matrix (RDM) of some spin block S is equal to the RDM in the 
state that time evolves with the mean-field Hamiltonian: 

ps{t) = tTs[e-^^^ |To) (Tol = tTs[U{t) l^o) (To| U\t)]. (C.28) 

Proof This is a direct consequence of Corollary C.7. 

Corollary C.9 (viz. Corollary 5.5) Let H E £ and |T) a state with cluster 
decomposition properties. // |T) is an excited state of the corresponding mean-field 
Hamiltonian 

Hl^\m) = E^\^) , (C.29) 

the expectation value of local observables in |T) is independent of time. Therefore, 
|\h) behaves locally as an excited state of H. 

The reverse is also true. If an excited state of H is locally equivalent to a state with 
cluster decomposition properties, then the latter is an excited state of the corresponding 
mean-field Hamiltonian. 

Proof Clearly the mean-field Hamiltonian is the solution of (5.3). Being 

|vl>) oc IT), by Corollary C.7 the expectation value of local observables is 
independent of time. The reverse holds true for analogous reasons. 

D. Self-consistency check of condition (4.8) 

Here we show that neglecting the anomalous term {cf (B.6)) in the time averaged 

Hamiltonian is a self-consistent approximation. To this aim, we consider the time 
evolution of a Slater determinant |To) under the Hamiltonian H, obtained from H 
by removing L.4,|“|. In Section 5 and Appendix C we proved that, as long as (9 is a 
local operator (but the class of allowed operators is in fact larger), can be replaced 
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by the mean-field time evolution operator U{t) (5.4). Here we show that inserting 
back at time t does not change the expectation value of local observables. In other words 
we are going to prove 

hm L{^o\UHt)[A\'^lO]U{t)\^o) = 0 Vt, (D.l) 

L—>-oo ^ J 

where (9 is a generic local operator. 

Using the notations of Appendix B, any local operator can be written as a linear 
combination of operators of the form 

^ = - (D.2) 

{pi} 

where 

n 

= n ■ (D.3) 

i=l 

If n is odd then (D.l) is trivially satisfied because U{t) |\['o) is a Slater determinant by 
assumption and hence the expectation value of an odd number of fermions vanishes. We 
therefore focus on the case n = 2m. The anomalous term of (B. 6 ) can be written 
as follows 

^ 1:^2 fc 3 ,/ C 4 

where and ^4 are the anomalous solutions of system (B.12), i.e. they are implicit 
functions of ki and k 2 defined by the system (B.12) and in addition fulfilling 


ki,2 ± ^ 3,4 7^ 0 , /C3 ± /C4 7 ^ 0, 


(D. 6 ) 


almost everywhere. Since |Tt) = U{t) ITq) is a Slater determinant, we can use the Wick 
theorem to compute expectation values. We then have 




iW 


w, 


iW. 


where contains terms in which two of the b’s in are contracted together 

and the other two are contracted with two b’s in O] C 4 [^|“|(!l]t contains all the terms 
in which any b in is contracted with a b in (9. 

According to the definition of any Wick contraction among b’s in it gives zero 
(because of (D.5)), hence the only non zero contribution to (D. 6 ) arises from C 4 [^|“|(! 2 ]t. 

To conclude the proof we will show that the terms in C 4 [^|“|C>]i scale as 0{L~^) and 
in the thermodynamic limit their contribution can thus be neglected. To this end it is 
sufficient to consider a typical element of C 4 [^|“j(!l]t 


;^;^G'rM])'^V(^i)b/32(^2)b/33(fc3)b/34(/c4)b^i(pi)b^2(P2)b<,3(p3) ■ ■ ■ b^3^(p2m), (D.7) 

where sums over the momenta {pi}, {ki} and the indices are understood, and we defined 
^a{p)hy{q) = {^t\ha{p)hg{q)\^t) (D. 8 ) 
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Ph}, 


(D,9) 


The 2m + 2 sums over the momenta are reduced to m by the Kronecker deltas arising 
from the Wick contractions. Because the number of factors L~^ exceeds by one the 
number of sums, the term turns out to be 0{L~^). The validity of equation (D.l) is 
then established. 


E. Additional properties of (7.1) 

E.l. Ground state and maximal energy state 

Here we show that the TFIC ground state iV'go)! with go satisfying (7.16), is the Slater 
determinant that minimises the energy (7.11) of the Hamiltonian (7.1). 

The symbol of the correlation matrix of the most general reflection symmetric one- 
site shift invariant Slater determinant is 

(E.l) 


r(/c) = nx{k)a^ -|- ny{k)a'^ -|- nz{k)a^ , 

with na{k) real functions that satisfy 

nl + nl + nl = l (E.2) 

nx,z{k) = -nx^z{-k) (E.3) 

ny{k) = Hyi-k). (E.4) 

Equation (E.2) manifests the fact that the initial state is pure, which indeed implies 
that the eigenvalues of r(A;) are ±1. The other conditions simply mean that the 
correlation matrix is a purely imaginary skew-symmetric matrix. The absence of a 
term proportional to the identity in (E.l) is a consequence of reflection symmetry. 

The energy (7.11) can be written as follows: 


£ = A 


^ny{k) 


TT 


9 


r\k (\k 

— ny{k) + / —[cos k ny{k) — sink nx{k)). (E.5) 


TT 


TT 


(E.6) 


The minimisation can be worked out by zeroing the variation of the functional 

n dk 

^nx,ny,nz,ij] = e - / — pik){nl{k) + nGk) + nl{k)) 

Jo tt 

with respect to its arguments. We immediately see that the variation with respect to 
Hz results in 

yi{k)nz{k) = d . (E.7) 

If we assume /i(A:) ^ 0, then nz{k) = 0 and we can enforce (E.2) by setting Ux = — sin^ 
and Uy = cos 6. Instead of working with the functional <I> we can therefore express the 
energy as 

/ dk \ 2 n dk n dk 

e[6'] = A( / —cos^fc) —g —cos6'fc-|- / —(cos/c cos 6^^-I-sin A; sin6'fc) (E.8) 

^Jo ^ 2 Jo TT Jo ^ 

and consider its variation with respect to dk- We obtain 

dp 


g — 2\[ — cos 6, 

0 


cos k 


sin 9k + sin k cos 9k = d 


(E.9) 
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which is solved by 


9o - 

^/l + g^- 2gQ cosk ’ 


(E.IO) 


where 5^0 satisfies (7.16). One of the solutions of (7.16) minimises the energy. The energy 
is instead maximal for 

= (E.ii) 

^Jl + gl- 2gi cosk 

where gi is a solution of 


9 = 9.-2a(^ 


dp gi — cosp 

’T yi+gf- 2t/i cosp 


(E.12) 


We emphasise that numerical data obtained by exact diagonalization of (7.1) for small 
chains are compatible with the two states of (E.IO) and (E.ll) being the actual ground 
state and the maximal energy state, respectively, of (7.1) in the thermodynamic limit. 


E.2. Excited states 


Following Corollary 5.5, we construct an infinite number of states for which the 
expectation value of local observables that time evolve under the Hamiltonian (7.1) 
is independent of time. To all intents and purposes, those states can be considered 
eigenstates of (7.1). 

A generic stationary solution of (7.7) should have y' = 0 and 

4{h — cos k)(j)k + 16(h cos k — l)yk = 0 , (E.13) 


namely 

- 16{h - cosfc)tr[r(fc)aV^'""] + 16(hcosfc - 1) tr[r(A;)a^] =0. (E.14) 


This is solved by 




—m{k) 


{hm — cos k)ay + sin k cr* 
-y/l + - 2hm cosk 


(E.15) 


where m{k) = m{—k) and we used ?/' = 0 to remove any term proportional to a^. The 
parameter hm is the solution of the equation 


hm 


g -\-2X 



dk hm — cos k 

—m{k )— , = . 

27r a/1 + - 2hm cosk 


(E.16) 


We recognise Tm{k) as the symbol of the correlation matrix of a reflection symmetric 
(Slater determinant) excited state of the TFIC model with magnetic field hm [67, 82, 
87,88]. Within this interpretation, if hm was independent of m, m{kj) = ±1 (for finite 
chains) would have provided an orthonormal basis of reflection symmetric excited states. 
We remind the reader that in the thermodynamic limit infinite excited states become 
locally equivalent and, in sufficiently regular cases, can be characterised by a function 
— 1 < m{k) < 1. However, in the latter characterisation, which exactly describes the 
local properties of the excited states, the manifest orthonormality properties of the basis 
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are lost. This problem is even more pronounced in our case, where depends on the 
excitation, and a naive hnite volume regularisation does not produce an orthonormal 
basis also for m{kj) = ±1 (this is not unexpected, being our description valid only in the 
thermodynamic limit). Nevertheless, the noninteracting states with correlation matrices 
(E.15) are locally equivalent to the excited states of (7.1). In addition, it seems that 
they span the entire space of reflection symmetric states. 

We notice that in the thermodynamic limit generally there are inhnite m that give 
rise to the same hm- Indeed we can always add to m{k) an even function Sm{k) such 
that — 1 < m + 6m < 1 and 



(E.17) 


which physically means the the excitation must have the same magnetisation per unit 
of length in the 2 ; direction. 

For the sake of completeness we report the expression of the energies of the excited 
states: 



(E.18) 


We conclude this appendix with a comparison between hnite chain data and 
expectations based on the mean-held correspondence, which strongly relies on the 
thermodynamic limit. In particular, we tried to apply Corollary C.9 to a hnite chain to 
estimate the lowest and highest excitations. Figures El and E2 show that the mean- 
held predictions are in fairly good agreement with numerical data even in small chains 
(notice that there can be 0{1/L) corrections in both directions). 
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